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FIGURES 


Comparison of isentropic ‘streamlines’ for; (a) bora (right-left flow), (b) chinook 
(left-right) and (c) Perth downslope winds (right-left). From Smith 1987, 
Brinkmann 1974, Pitts and Lyons 1989, respectively. 


Composited upwind and downwind soundings for chinook windstorm 
events (after Brinkmann 1974). 
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Map of Mount Lofty Ranges. 
Map of Adelaide and surrounds. 


Steady-state streamlines for stratified flow over a ‘Witch of Agnesi’ obstacle with 
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H = 500 m, a = 3 km, No = 0.01 s"!, and up = 5 m/s (from Laprise and Peltier 1989a): 


(a) a linear solution, (b) a nonlinear solution, (c) a nonlinear and nonhydrostatic 
solution and (d) a nonlinear and nonhydrostatic solution with a = 1 km. 


Analogue of hydraulic flow as suggested by Kiittner. 
Schematic diagram of the splitting-streamline concept of S85. 


Mean wind at Highbury (near Wattle Park) in January and February 1980 ona 
diurnal basis for west to southwesterly streams and for east to southeasterly 
streams (after Marsh 1987). 


Mean sea level pressure chart for 0900 CSST 13 February 1980: typical for 
strong gully-winds. 


Surface synoptic chart for 0900 CSST 7 February 1985. 


Synoptic charts for 0900 CSST 7 February 1985 for; (a) gradient, 
(b) 850 hPa, (c) 7000 feet, (d) 700 hPa, (e) 14 000 feet, (f) 500 hPa. 


Temperature profiles for Adelaide Airport and Mount Gambier at 0900 CSST on 
7 February 1985. ML denotes Mount Lofty summit. 


Charts for 0900 CSST 26 February 1991: (a) gradient level and (b) mean sea level 
pressure chart. 


Temperature profiles for Adelaide Airport and Mount Gambier at 0900 CSST on 
26 February 1991. When used in conjunction with Fig. 14, these show that there 
is a shallow layer of cool air crossing the ranges from the east to southeast. 


Mean sea level pressure chart for 2100 CSST 6 February 1972. 


Acoustic-sounding record of turbulent-layer depth, the vertical velocity at 100 m 
and the 100 m-layer mean wind for the period 2200 CSST on 6 February to 0200 
CSST 7 February 1972 at Salisbury, South Australia. 
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Vertical temperature profile for 2100 CSST 6 February 1972 at Adelaide Airport. 


Mean sea level pressure charts for 0900 CSST on; (a) 3 February 1988 and 
(b) 4 February 1988. There is weak gradient flow for (a) but a shallow cool-air 
intrusion occurs in the wake of the (dotted) trough in (b). 


Three-hourly isentropic cross-sections (using virtual potential temperatures 
in °C) for night of 2/3 February 1988 along the true bearing 110°-290° 
through Mount Lofty. Faint lines show aircraft tracks and balloon trajectories. 
Offset winds are for Adelaide Airport (left) and Nairne (right). 


As for Fig. 20 but for specific humidity in gm/kgm. 
As for Fig. 20 but for the night of 3/4 February 1988. 
As for Fig. 22 but for specific humidity in gm/kgm. 


Three-hourly temperature profiles at Nairne (solid) and Adelaide Airport (broken) 
on the nights of the 2/3 and 3/4 February 1988 (identified as nights A and B, 
respectively). At 2330 CSST on the second night a sounding was completed at 
the crest of the ranges. The superadiabatic layers on the B3 profile have been 
double checked and therefore appear to be real. 


Average AT as a function of Adelaide Airport wind direction at 736 m, above 
MSL (after Marsh 1987). 


Scatterplot matrix of peak gust at Port Stanvac (between 1000 and 1200 CSST), 
AT and geostrophic wind speed at 1100 CSST. 


Photographs (looking west out over Backstairs Passage) of hydraulic jump in fog 
bank on 2 May 1988 in Backstairs Passage, South Australia. Top sequence taken at 
1420 CSST; lower one taken at about 1530. See Grace (1991) for further details. 
(Photographs courtesy of Mr Richard Geytenbeek) 


Steady 2D flow over a mountain with upstream static stability No, and velocity, up. 
The streamline initially at altitude dp splits, so that above dy the flow is unchanged. 
Below the lower split streamline the flow is essentially hydraulic, and between 

the split streamlines the air is well-mixed and turbulent but with no mean motion. 
Stability and speed of the layer beneath the lower split streamline (69) are constant 
with height. 


Graphical solution of Eqn 18 showing do and higher order solutions as a function 
of A. Also shown is the ratio dH for n= 1. 


Plots of do, di, dy, dp, Ave and dj, against A. 


Plots of jump strength B, wind-speed amplification 3 and turbulence intensity 7 
against A. 
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Nondimensional drag, 7, for S85 (labelled S), Q1'/2 model (labelled Q) and 
linear waves compared as functions of A with spot values from numerical 
experiments by Durran and Klemp (1987) (triangles) and Bacmeister and 
Pierrehumbert (1988) (circles). 


Model solutions of do for Q1'/ and S85/SS87 as a function of fr with spot 
values from tank experiments (unfilled circles) and numerical experiments (filled). 
The spot values and S85/SS87 curve are from Rottman and Smith (1989). 


The defining sketch for M. The topography has been modified in such a way that 
the lee plains are lower than the windward ones. 


Plots of (a) j jump strength, B, (b) wind-speed amplification $ and (c) turbulence 
intensity j against A for various values of M. 


Plan view of a string of conical hills ‘melded’ together. Heavy solid lines are 
contours and dashed lines represent streamlines. Flow along VV represents flow 
through a valley while SS represents flow along the highest path or spur line. 


Plots of (a) jump strength, B, (b) wind-speed amplification 3 and (c) turbulence 
intensity j against A for various values of v 


Adaptation of the diagram of Houghton and Kasahara (1968) for transitional and 
nontransitional flow regimes to the Q1'/, model. 


Comparison of the results of Durran and Klemp (1987) for high drag (filled circles) 
and low drag (hollow circles) events with Q1'/ predictions. 


Plot of dual periods in minutes against H/ug and No. 


Figures 41 to 49 — evolution of Type I flow with ‘Witch of Agnesi’ topography 
given by H= 500 m,a=3 km, b= - km, and preconditions of up = 10 m/s, 
Vo = 0 mis, dy = 1500 m, g’ = 0.2 mls*, f =-83 x 10 ©s"!, C, = 0.002 and 
E=0.0 unless noted otherwise. 


The case f = 0 s"!, Cp = 0, in which transitional flow and downslope winds 
eventuate. 


The case in which Cp = 0 but f = -83 x 10° 5"!. In this, the leeward jump is 
further leeward but the upstream jump is weakened and retarded. 


The ‘reference’ preconditions for Type I flow. 
Non-transitional flow for the reference conditions except that u = 5 m/s. 
Non-transitional flow for the reference conditions except that dy) = 2000 m. 


Flow results for the reference conditions except that g’ = 0.1 m/s. Note that the 
supercritical flow is more extensive than before, although maximum speeds are 
decreased. 


Page 


40 


40 


42 


43 


44 


46 


47 


47 
52 


59 


60 
60 
61 
61 


63 


Figures (continued) 


Figure 


47 


48 


49 


50 


51 
52 
Al 
Bl 
Cl 


Page 


Flow results for the reference conditions except that g’ = 0.3 m/s”. Note that the 
transitional flow is much less extensive, although maximum speeds are increased. 
This is equivalent to very strong gully-winds about the foothills with much weaker 
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Flow resuits for the reference conditions except that Cp over the hills region is 
increased to 0.01 (though Cp over the plains is held at 0.002). Note that the 
leeward jump is held back. 
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leeward jump is held back, but less so than in Fig. 48. 


Type II flow with initial conditions of a southerly of 8 m/s, dy = 400 m, g’ = 0.1 mls? 
and Cp = 0.002. Note that this type of flow is associated with a sea-breeze from the 
Coorong area. 


Type II flow under the same conditions as for Fig. 50, except that vo = -4 m/s. 
Type II flow under the same conditions as for Fig. 51, except that E = 0.001. 
Schematic diagram of transitional flow over an obstacle. 
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Schematic diagram of homogeneous fluid crossing a two-dimensional obstacle. 
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SYMBOLS AND ABBREVIATIONS 


> 


denotes non-dimensionalised value 
leeward position where h = 0 


jump strength 


<~ WR 


measure of valley constriction 


AT temperature difference 


€ amount of perturbation 

6 potential temperature 

K numerical diffusion coefficient 

Xr numerical space-time grid ratio 
et ML times U = udu; a Momentum flux 
p density 

o N times U = vdu; a momentum flux 


a,b mountain half-widths 

BP — Bacmeister and Pierrehumbert (1988) 
c horizontal gravity-wave speed 

Cp _ coefficient of surface friction 


CST Central Standard Time 
(= UTC + 9'/) hours) 

CSST Central Standard Summer Time 
(= UTC + 10'/: hours) 

d layer depth 

DK _ Durran and Klemp (1987) 


e energy loss at jump 

E coefficient of interfacial drag 
F Froude number 

f frequency, Coriolis force 

g gravity 

g modified gravity 

h topographic height 

H mountain height (maximum) 
HK — Houghton and Kasahara (1968) 
J turbulence intensity 

m cross-range flux (ud) 


M height ditference of plains compared to H 


MF — momentum flux 


n along range flux (vd), dummy index 

N _ Brunt-Vaisala frequency 

P pressure 

P period 

PC Peltier and Clark (1977) and subsequent 
papers 

PF _ pressure force 

m mass flux, specific humidity 

Ql'/2 steady-state quasi 1'/-layer model 

QT time-dependent version of Q1'/: model 

r mountain-wave drag 

$ ratio of windward to leeward wind 

S85 Smith (1985) 

SS87 Smith and Sun (1987) 

T temperature 

u cross-range wind speed 

v along-range wind speed 

l Niu 

l, Scorer parameter 

Subscripts 

0 upstream, upper 

1 before jump 

Ps after jump 

A aloft 

é crest 

S S85, or Scorer 

1 i threshold or transition 


ABSTRACT 


Downslope winds occur in the lee of all mountain ranges. Those occur- 
ring over the western flank of the Mount Lofty Ranges of South Austral- 
ia are known locally, but slightly misleadingly, as gully-winds. Case and 
field analyses of several gully-wind events are presented. Within the 
framework of hydraulic theory, downslope winds and, in particular, the 
Adelaide gully-wind are investigated and modelled. A steady-state mod- 
el for continuously stratified flow perturbed by an obstacle is developed 
and compared with other models. Comparison with the results of numer- 
ical model and tank experiments shows good agreement. A useful result 
of this steady-state model is that it is apparently reasonable to approxi- 
mate a continuously stratified layer by an equivalent homogeneous layer. 


With this approximation, and with some findings from the field and case 
events as to the scale depth of the gully-winds, the time-dependent shal- 
low-water equations are used to model two proposed archetypes of Ad- 
elaide gully-wind situation: a Type I, where a synoptic stream is 
perturbed by topography; Type II, where a density current flowing ini- 
tially parallel to a mountain range is subject to inertial forcing. A sur- 
prising feature of Type II is that Coriolis turning is sufficient to force 
the density current to spill over the range as a supercritical leeslope 
flow. Surface drag is found to be important. Results of Type I and II sim- 
ulations compare favourably with case studies compiled from synoptic 
data, with field studies and with sophisticated numerical models. 


Using the results of the steady-state model, a resonance mechanism for 
gustiness is proposed and developed. The susceptibility of supercritical 
(single-layer) flow to perturbations is investigated and it is shown that 
supercritical flow would tend to jump near the lee foot. It is concluded 
the relatively simple hydraulic theory developed herein can explain 
most of the features of the Adelaide gully-wind. 


The time-dependent model has been developed for operational use in a 
forecasting office. Rigorous and objective testing is precluded because 
of inadequate upstream observations and continuously changing synop- 
tic features. Nevertheless, predictions tend to be consistent with fore- 
casters’ subjective assessments. The software package is included with 
the thesis. Three published papers (one co-authored) are also presented 
as part of this thesis. 
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CHAPTER 1 


INTRODUCTION 


1.1 General 


Downslope winds occur in the lee of mountains all over the world. Well-known downslope winds 
include the Féhn of Alpine Europe, the Mistral of the northern Mediterranean coast, the bora of 
Yugoslavia, the chinook of North America and the zonda of Argentina. Less widely known are 
the Norwegian fallwind (Humphreys 1929), the helm wind of the British Pennines (Manley 
1945), the strong ‘katabatics’ off the coastal iceshelves of Antarctica and Greenland (Ball 1956) 
and those down the eastern flank of the Western Cordillera of Colombia (Lopez and Howell 
1967), the hiroto-kaze of Japan (Suzuki and Yabuki 1956), the chinook arch of the Canadian 
Rockies (Ash 1987), the scarp wind of Perth (Van Burgel 1986; Pitts and Lyons 1989,1990) and 
the Adelaide gully-wind (Bureau of Meteorology 1971; Holton 1985; Marsh 1987; Grace and 
Holton 1988,1990). One of the most spectacular and most studied downslope winds is that of 
Boulder, Colorado (e.g. Colson 1954; Lilly and Zipser 1972; Brinkman 1974; Klemp and Lilly 
1975; Peltier and Clark 1979; Bower and Durran 1986; Durran 1986). 


In some locations, downslope winds may gust to more than 50 m/s, creating a hazard to 
aviation in particular and damage to buildings and trees. The bora has even overturned streetcars 
in Trieste! In one famous incident, the hangar of the research aircraft, that was on a sortie studying 
the chinook, was severely damaged (Lilly and Zipser 1972)! Typical, but not universal, features 
of such winds are a night-time speed enhancement, sudden onset, a threshold speed to windward 
(Nicholls 1973; Brown 1986), marked variability in strength but not in direction, and a stable lay- 
er or inversion near or just above hilltop level (Atkinson 1981). Generally the cross section of the 
mountain range is ramp shaped — that is, gently sloping to windward with a much steeper slope to 
leeward. 


Much investigative effort has gone into compiling a quantitative description of downslope 
winds. Although the bora is famous as a cold wind and the chinook as a relatively warm wind, the 
mechanism for each is probably the same and the relative warmth of the chinook is due to the dis- 
placement of a very cold shallow layer over the lee plains. There are also many counter-examples 
of warm bora and cold chinooks. Boras and chinooks have been known to contain gusts as high 
as 50-60 m/s. Downslope winds are usually more pronounced at night. Typically they are very 
gusty yet very steady in direction. For climatologies of chinooks and boras see Brinkmann (1974) 
and Yoshino (1976). 


Isentropic cross-sections from a bora (Smith 1987), a chinook (Brinkmann 1974) and 
Perth’s downslope wind (Pitts and Lyons 1989) are compared in Fig. 1. (In passing, note that 
Smith (1987) provides a series of magnificent colour photographs of bora effects). All three, and 
many others in the literature, have in common a single large amplitude (forced) wave with or 
without secondary lee waves further downstream. Turbulence and vertical motion are most pro- 
nounced in association with the large trough. Most of the ‘action’ is within the lowest 2-3 km, 
although this is not always so: Lilly and Zipser’s (1972) famous example, featured in many text 
books, shows that the whole tropopause descended to within about 2 km of the leeside surface. 


Brinkmann (1974) composited upwind and downwind soundings associated with Boulder 
chinook windstorms and the results are shown in Fig. 2. Referring to the latter, it can be seen that 
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Fig. | Comparison of isentropic ‘streamlines’ for; (a) bora (right-left flow), (b) chinook (left- 
right) and (c) Perth downslope winds (right-left). From Smith 1987, Brinkmann 1974, 
Pitts and Lyons 1989, respectively. 
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the inversion is consistently lower to leeward. For the bora, it is known that a windward inversion 
occurs somewhere between 1500 and 2000 m (AMSL), whereas to leeward the inversion is usu- 
ally some 500 to 100 m lower (Yoshimura 1976). This is very similar to Suzuki and Yabuki’s 
(1956) findings for a hiroto-kaze (a downslope wind event accentuated by valley effects). They 
presented a transect of soundings for five stations, showing a pronounced inversion decreasing 
from an altitude of 3000 m on the windward side to 1600 m on the leeward side of the ranges. 
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Fig. 2 Composited upwind and downwind soundings for chinook 
windstorm events (after Brinkmann 1974). 


1.2 The Adelaide gully-wind 


The Mount Lofty Ranges of South Australia have a ramp shape with respect to easterly flow. 
These ranges extend northeastward from Cape Jervis to Myponga, then approximately north- 
wards on past Adelaide. The ranges are 40 to 50 km wide with the highest point being Mount 
Lofty itself which is 736 m AMSL. The average height of the ridge is about 500 m, however. To 
the east, the countryside is generally flat and open and about 100 m AMSL; to the west, the land 
falls steeply over a horizontal distance of 3 to 5 km, but then slopes gently across a coastal plain 
of varying width. (Figures 3, 4 and 5 show locations and topography.) Further north, the northern 
Flinders Ranges have a ramp shape to westerly flow. 


Easterly downslope winds commonly occur on the western flank of the Mount Lofty Rang- 
es and often extend several kilometres further westward. It is possible that, on occasion, they blow 
right across Gulf St Vincent although this conjecture is based on undocumented evidence from 
local fishermen. In the Adelaide region, these winds are known as gully-winds and a common ob- 
servation is of moderate easterly winds at the mountain top with strong gusty easterly winds about 
the lee foothills, while a few kilometres further downstream the wind is light and may even be 
weakly reversed. Frequently the wind above 1500-2000 m on a synoptic scale is westerly. In ex- 
treme cases trees have been uprooted and houses damaged. 
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Fig. 3 Map of South Australia. 


The winds may seriously hamper fire-fighting operations in the Mount Lofty Ranges. 
Marsh (1987) reported a gust of 30 m/s at a wind energy monitoring site at Myponga Beach 
which is located about 40 km south of Adelaide. Holton (1985) documented aircraft reports of 
turbulence at Adelaide Airport and found that 74 per cent of all reports rated as severe occurred 
during gully-wind events. Major errors (> 5°C) in summertime minimum temperature forecasts 
for Adelaide are mainly attributed to the unpredictable nature of the gully-winds (Furler private 
communication). Late last century, during the era of Adelaide’s horse-drawn trams, extra horses 
were required for trams travelling from Parkside (a suburb on the southern edge of the Adelaide 
city centre) to Glen Osmond into the face of strong gully-winds (Kingsborough 1967). Typically, 
they occur if the synoptic stream is a stable and vigorous east to southeasterly. These conditions 
correspond well with the pre-conditions generally accepted for forecasting large-amplitude 
mountain waves. Durran (1986a) lists these as being stability at lower levels, wind speed of at 
least 7 m/s (for mountains of | km) and asymmetric mountains with gentle windward and steep 
leeward slopes. 


The Adelaide region downslope winds are more pronounced in gully sites, hence the local 
term ‘gully-wind’. On occasion, despite apparently favourable pre-conditions, downslope winds 
in the lee of the highest part of the ranges around Mount Lofty are weak or non-existent, yet 10 to 
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Fig. 4 Map of Mount Lofty Ranges. 


50 km north or south, where the ranges are about 200 m lower, strong downslope winds are ob- 
served. Little is documented of downslope winds elsewhere in South Australia: during an envi- 
ronmental impact study (Holmes 1981) in the northern Flinders Ranges, westerly downslope 
winds of 28 m/s (mean) were recorded: this compares to 30 m/s for the hangar-damaging chinook 


incident. 


1.3 Philosophy of approach 


In keeping with Occam’s Razor, this thesis pursues the simplest and least abstract (in the author’ s 
opinion) of the postulated mechanisms of downslope winds, and extends it to include more real- 
istic model atmospheres. Extremely simple models have the advantage that if the ‘correct’ be- 
haviour is observed or obtained then it is likely that the underlying physics has been captured. 
Most mesoscale models include elaborate parametrisations of physical processes (that is, bound- 
ary-layer formulation, boundary conditions, humidity and drag) and highly complex model struc- 
ture; it is sometimes difficult to distinguish between physical and non-physical aspects of the 
model results (Hoinka 1985). The particular theory applied here has much in common with 
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observations of water flow both in the laboratory and in nature and is therefore easily visualised. 
An associated aim of this work is to produce a workable model on a modest desktop computer for 
real-time use by operational forecasters. 


1.4 Outline of paper 


Chapter 2 summarises current models and theories of downslope winds. Chapter 3 presents case 
and field studies of the Adelaide gully-wind, and shows photographic evidence of a hydraulic 
jump in a bank of sea fog at the southern end of the Mount Lofty Ranges. Chapter 4 develops a 
steady-state ‘layer-and-a-half’ model and compares its results with those of other models in the 
literature. A possible cause for the quasi-periodic gustiness of downslope winds is presented in 
Chapter 5. Chapter 6 develops a time-dependent model of downslope winds. Results are com- 
pared to available observations. The importance of drag is examined. The discussion of Chapter 7 
draws together the literature, the data presented, and the model results. Chapter 8 points to some 
avenues for further research. Figures 1 to 3 show relevant topography and all South Australian 
locations as mentioned in the text. Because the Adelaide gully-wind is a markedly diurnal sum- 
mertime feature it is natural to use CSST (Central Standard Summer Time). CSST = UTC+10", 
hours. On one or two occasions Central Standard Time is used, where CST = UTC+91/, hours. 


CHAPTER 2 


THEORY AND CURRENT MODELS 


Although the mechanism for downslope winds is still not established (Scorer 1978; Atkinson 
1981 and Durran 1986 a,b), it is well-established that their ultimate cause is the buoyancy force 
in conjunction with topography. Currently several different mechanisms are invoked to explain 
the development of downslope winds. The review contributions of Durran (ibid.) are particularly 
lucid and comprehensive. Corroboration of each of the theories is mainly attempted through the 
use of numerical models incorporating the fundamental equations and these in turn have been val- 
idated against a few well-documented case studies, or the analytic solution of idealised atmos- 
pheres, or laboratory tank models. The tank model experiments have also been used to test 
implications of the postulated mechanisms, as have the available observational data. 


2.1 Background: wave theory 


The perturbation of incident hydrostatic flow by the gravity waves launched by a ridge was inves- 
tigated by Queney (1948). Two types of wave are induced by a ridge — a forced wave which is 
collocated with the topography and lee waves which propagate downstream (Pielke 1984). 
Trapped lee waves are, of course, acommon sight. Scorer (1949) showed that the lee waves would 
be well developed if the Scorer parameter, /,, decreased with height where 


2 
N 1924 
l= |>--s qd) 
ue uae 


where N is the Brunt-Vaisala frequency, u is the cross-ridge component of the wind speed and z 
is altitude. Nevertheless, observational experience (e.g. Vergeiner 1971; Brinkmann 1974; Pitts 
and Lyons 1989) and numerical simulations (e.g. Hoinka and Rosler 1987) show that the trapped 
lee waves are not directly relevant to downslope winds. Downslope winds are associated with an 
amplification of the forced wave. Most of the numerical models used have been hydrostatic but 
this is not a serious limitation, as they do simulate the forced wave quite well, although they do 
not properly simulate the trapped lee waves (Pielke 1984). 


For a nonhydrostatic regime, Long (1953) showed that if the upstream profiles of stability 
and speed (No and up) are constant with height then the streamline deflection 6 satisfies the linear 
equation (Long’s equation) 


vs 08 a 
eet ald = 0 (2) 


where x is the direction of stream flow and / is a specialised form of Scorer’s parameter; namely, 


Pe ae 
Ug (3) 


If the hydrostatic approximation is made — which is equivalent to assuming that 
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This linear perturbation equation has a solution 


5 = he’ (6) 


where hi is a complex function whose real part defines the surface topography, h(x). (That is, 


Rth = h). For a bell-shaped (‘Witch of Agnesi’)! mountain of the form 


H 





=| 
Il 


3 7) 


with H and a being mountain height and ‘half-width’ respectively, then the solution is given by 


5 = 





1x (8) 
an 


The form of Eqn 7 was originally chosen by Queney because it is simple, realistic and has a con- 
venient Fourier integral of the form 


= Hal e™ cos (kx) dk (9) 
0 


where k is wave number. The use of the ‘Witch of Agnesi’ profile has since become traditional in 
downslope wind models. 


It turns out that for a symmetrical Witch-shaped mountain critical steepening of the stream- 
lines occurs if /H = / (see Fig. 6(a)). By critical steepening, we mean that streamlines become 
vertical, on the verge of overturning. The lower boundary condition used is that § = A on z = 0 
(not on z = A, the actual surface) and this effectively assumes that the mountain is a ‘small’ per- 
turbation”. This leads to a linear theory, in that the amplitude of the forced wave is linear, and 
hence the term ‘linear’ waves. The lower boundary condition is clearly not realistic and results in 
the lowermost streamlines going through the mountain. Note that the dashed line in Fig. 6(a) rep- 
resents the ground surface. 


1. In the theory of radiation, this function is sometimes known as the Lorentz function, so it could be described 
as a ‘Lorentz-shape'’. 


2. This condition is sometimes known as the slender-body approximation in aerodynamics. 
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Fig. 6 Steady-state streamlines for stratified flow over a ‘Witch of Agnesi’ obsta- 
cle with H = 500 m, a = 3 km, No = 0.01 |, and uy = 5 m/s (from 
Laprise and Peltier 1989a): (a) a linear solution, (b) a nonlinear solution, 
(c) a nonlinear and nonhydrostatic solution and (d) a nonlinear and non- 
hydrostatic solution but a = 1 km. 


{ non-linear theory requires Eq. 5 tc be solved subject to the non-linear boundary con- - 
dition that 0 = h at z = h. For detailed discussions of the linear and nonlinear theories, the reader 
is referred to Gill (1982) or Laprise and Peltier (1989a). In the nonlinear case, critical steepening 
occurs if IH = 0.85 (see Fig. 6(b)). Baines and Granek (1990) have shown that critical steepening 
occurs for 0.5 < /H < 1 where the lower (higher) value applies to ramp- (backward ramp-) shaped 
‘Witch of Agnesi’ mountains. 


In the linear theory, critical steepening first occurs at an altitude of (n+3/4) of a vertical 
wavelength, where n = 0,1,2..., and is directly above the crest. An implication of having the max- 
imum steepening over the crest is that the strongest winds would occur halfway along the lee 
slope; whereas, in fact, they usually occur near the foot. 
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In the nonlinear theory, the altitude of maximum steepening remains at about the same height 
and is displaced downstream (with some dependence on the half-width of the mountain). Laprise 
and Peltier (1989a) have investigated the nonlinear nonhydrostatic case and the essential features 
of their solutions compare closely (Fig. 6(c)). For narrower hills, the hydrostatic approximation 
is coarser and secondary lee waves become more prominent (Fig. 6(d)). 


2.2 Laboratory models 


These are usually tank models of saline water in which the density decreases linearly with height. 
A model mountain is towed from a standing start and a camera travelling alongside records the 
fluid motion made visible by dye streaks or neutrally buoyant polystyrene balls. Experiments with 
tank models (e.g. Baines and Hoinka 1985; Boyer and Tao 1987; Rottman and Smith 1989) have 
shown conclusively that the flow is dominated by lee waves for small H/ but for larger H/ the flow 
evolves to one of hydraulic character. 


The advantages of the tank models are that they show the extreme complexity of the phe- 
nomena and they can be used as a check on numerical models directly and to test implications of 
the mechanisms postulated. A disadvantage is that some effects such as boundary-layer drag may 
be disproportionately more or less important in the atmosphere. Currently, there are three mech- 
anisms postulated to explain downslope winds. 


2.3 Hydraulic mechanism 


The first mechanism is often referred to as that of the hydraulic jump (Kiittner 1939, 1958; Long 
1953, 1954; Ball 1956; Lopez and Howell 1967; Houghton and Kasahara 1968; Arakawa 1968), 
where an inversion (usually) above ridgetop height is regarded as analogous to the upper surface 
of a shallow layer of water. 


The hydraulic jump model of downslope winds is based on the premise of a transition from 
subcritical windward flow to supercritical leeward flow which then dissipates its energy within a 
hydraulic jump. On both windward and leeward sides, the flow accelerates (prior to the jump). 
Kiittner’s idea is represented in Fig. 7. This approach has the merit of explicitly accounting for 
the nonlinearities in the system; however, unlike the atmosphere, all the models considered have 
an upper boundary consisting of either a free surface or a rigid lid. In this type of model, no up- 
ward energy transport can occur across the upper boundary, and so vertically propagating atmos- 
pheric waves cannot be properly represented. Because many of the hydraulic models have been 
limited to homogeneous layers, any comparison of results with waves in a (more realistic) contin- 
uously stratified atmosphere must be primarily qualitative. This thesis seeks to explain much of 
the associated observations by applying and extending the hydraulic theory. 


2.4 Lee wave amplification — Klemp and Lilly 


A second mechanism for downslope winds is based on an extension of Scorer’s linear lee wave 
theory by Klemp and Lilly (1975, 1978), who proposed that amplification of lee waves (primarily 
the forced wave) could be accomplished by the resonant partial reflections of wave energy prop- 
agating vertically between the surface and an upper reflecting level, or levels. This mechanism is 
sometimes referred to as linear lee wave amplification. 
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Fig. 7 Analogue of hydraulic flow as suggested by Kiittner. 


The upper reflection occurs at a discontinuity in /. Strictly speaking, the reflection occurs at 
discontinuity in Scorer’s parameter, /,, where /, is given by 


2. 2 
5 
peg ee Dd (10) 


Such a discontinuity is often termed an Environmental Critical-Layer (ECL) which is usually 
modelled by a wind-shear layer one or two kilometres deep; less commonly, a strong inversion 
layer is used. The implication is that if the upstream atmosphere is ‘tuned’ in such a way that, for 
example, the tropopause is located at an altitude equal to a half wavelength (or some integral mul- 
tiple thereof) of the internal (vertical) waves, then downslope winds will occur. 


However, since this is a linear model, it must be applied with caution to the case of large- 
amplitude waves, because the perturbation fields associated with these waves are comparable in 
magnitude to the mean flow. Predictions from linear multilayer models are sometimes rather mis- 
leading ((Durran 1986a,b); and Durran and Klemp (1987) and Bacmeister and Pierrehumbert 
(1988), (hereafter DK and BP respectively)) found little support for this mechanism. 


2.5 Lee wave amplification — Peltier and Clark 


A third mechanism for the generation of large-amplitude waves has been suggested by Peltier and 
Clark (Clark and Peltier 1977, 1984; Peltier and Clark 1979; Laprise and Peltier 1989a,b,c). They 
examined the behaviour of large-amplitude mountain waves with a numerical nonlinear non-hy- 
drostatic model which gave a detailed representation of the atmospheric flow. They observed that 
a substantial increase in the strength of the wave occurs once the streamlines overturn and the 
wave ‘breaks’. The wave-breaking region is characterised by strong mixing and a local reversal 
of the horizontal wind. They postulated that the energy in the upward-propagating wave is trapped 
below this self-induced critical-layer (the region where the wind reverses direction), thereby pro- 
ducing large-amplitude waves. This critical-layer developed if the local Richardson number de- 
creased below 0.25. 
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Fig. 8 Schematic diagram of the splitting-streamline concept of S85. 


At present, there is still some controversy about the role played by self-induced critical-lay- 
ers on the subsequent amplification of the wave. An atmosphere ‘primed’ with an environmental 
critical-layer at an altitude equal to (n+3/4) of the vertical wavelength should. by this theory, pro- 
duce downslope winds. Though some numerical and tank experiments support this prediction, 
many others do not. Smith (1987) has suggested that the data from such experiments support an 
internal hydraulic mechanism, rather than a mechanism based on vertically propagating waves. 
However. more recently, Laprise and Peltier (1989a,b,c) have concluded that the ‘linear instabil- 
ity of the nonlinear wave field’ explains the transition of the flow to one with large amplitude 
waves, 


2.6 Smith’s model 


In a landmark paper, Smith (1985), (hereafter S85), and subsequently Smith and Sun (1987), 
(hereafter SS87), noted that the atmosphere above large ampiitude waves remains relatively un- 
disturbed and used this idea to extend Long’s model to continuously stratified fluid. On the as- 
sumption that there exists a decoupling of the low-level accelerated flow (termed the 
bydraulically active Jayer by him) S85 developed a steady-state model in which the wave energy 
is trapped within the low-level flow (Fig. 8). The height of the initial decoupling, /dp, is depend- 
ent on the non-dimensional mountain height JH. 


Numerical experiments by DK and BP showed that in the presence of an environmental critical- 
layer near or below Smith’s do, lee-wave amplification occurred, and the overall flow was remarka- 
bly similar to hydraulic flow. For an environmental critical-layer located above dp, they found that 
similar transitional flow could occur provided /H > 0.85 — which is the well established condition of 
Miles and Huppert (1969) for wave breaking in an unbounded medium. However, even with /H < 
0.85 and with an environmental critical-layer located well above do, but close to dy + 2nz7/, transi- 
tional flow could occur, although the resonance was considerably weaker and required ‘tuning’. 
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BP found that the character of the high-drag states compared well with the hydraulic analysis 
of S85 and could not be explained in terms of linear resonance (Klemp and Lilly’s postulated mech- 
anism); and further concluded that regardless of the appropriate upper-boundary condition (that is, 
whether the environmental critical-layer is reflective or absorptive in nature), it is clear that nonlin- 
ear effects captured by hydraulic theory are the key element. Only partial support for the validity of 
PC’s mechanism was found; for example, flow with an environmental critical-layer at the resonance 
altitude of Clark and Peltier (1984) did not necessarily develop to transitional flow as expected. Pitts 
and Lyons (1990) found partial support for PC’s proposed mechanism by comparing numerical mod- 
el results with their observations of Perth’s downslope winds. Consistent with BP, DK concluded 
that their results suggested that ‘the strong lee slope winds associated with wave-overturning are 
caused by a continuously stratified analog to the transition from subcritical to super-critical flow in 
conventional hydraulic theory’ and, further, that their results did not appear to be consistent with the 
amplification mechanism of PC. Rottman and Smith (1989) state explicitly that only the model of 
S85 (as extended by SS87) is consistent with all published results of numerical model and tank ex- 
periments. 


2.7 The model of Baines and Granek 


S85 assumes pre-existing decoupling between hydraulically active and passive layers and that the 
decoupling is due to either the wave-breaking effect of PC or, for subcritical mountains, an environ- 
mental critical-layer positioned ‘near’dy + 2n7/I (n = 0.1,2,...) although the exact meaning of ‘near’ 
is not specified. Baines and Granek (1990) have extended this concept and formulated a unifying, 
comprehensive and quantitative steady-state model. 


They provide nonlinear wave-field solutions of Long’s equation for asymmetric ‘Witch of Ag- 
nesi’ mountains up to the nondimensional height of 47. which they show to be a function of asym- 
metry. (For convenience, the nondimensionalised height is denoted as Al and defined as /H; and Hr 
denotes the threshold value of A at which streamline ov erturning first occurs.) For H = Ar. the so- 
lution is that of S85. They assume that asf is increased beyond Hy, the upper surface streamline 
remains unchanged: this assumption being supported by the experiments of Baines and Hoinka 

(1985). For the region H > Ay. a multilayer (64 level) steady-state hydraulic model based upon the 
numerical solution of the shallow-water equations is used. Beyond about H ~ 1.5. blocking of low- 
er-level fluid occurs and the effective value of A is changed. This overall approach accords with sug- 
gestions by Smith (1987) and Rottman and Smith (1989). 


CHAPTER 3 
OBSERVATIONS OF ADELAIDE GULLY-WINDS 


3.1 General 


Adelaide gully-wind events occur mostly in the warmer months of November to March and dur- 
ing the night (Bureau of Meteorology 1971). Marsh (1987) analysed 8 days of west to southwest- 
erly airstream and 18 days of east to southeasterly airstream for January and February of 1980. 
Using data from a Woelfle anemometer at a quarry site in Highbury (where an Environmental Im- 
pact Study was conducted), mean diurnal winds for each stream type were compiled (Fig. 9) and 
the difference between them was found to be quite remarkable: at Highbury southwesterly 
streams decayed nocturnally, as one might expect, whereas east-southeasterly streams freshened 
nocturnally to 8 to 9 m/s. 


Marsh (1987) identified the most common synoptic situation for gully-wind events: namely, 
a high centred at 40°S between 125° and 140°E, with a cold front moving across or east of Tas- 
mania and extending a weak northern portion into the southeast of South Australia. This north- 
ernmost portion tends to become quasistationary and decay. An example is given at Fig. 10, which 
shows the synoptic situation at 0900 CSST on 13 February 1980. Of course, local forecasters have 
been aware of the significance of these patterns for many years. 
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Fig. 9 Mean wind at Highbury (near Wattle Park) in January and February 
1980 on a diurnal basis for west to southwesterly streams and for 
east to southeasterly streams (after Marsh 1987). 
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Fig. 10 Mean sea level pressure chart for 0900 CSST 13 February 1980: 
typical for strong gully-winds. 
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Fig. 11 | Surface synoptic chart for 0900 CSST on 7 February 1985. 
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Fig. 12 | Synoptic charts for 0900 CSST on 7 February 1985 for; (a) gradient, (b) 850 
hPa, (c) 7000 feet, (d) 700 hPa, (e) 14 000 feet and (f) 500 hPa. 


Grace and Holton (1988) examined 52 easterly events (geostrophic flow from the northeast 
to southeast quadrant at night time over Adelaide) occurring between 1986 and 1987. Thirty of 
these were gully-wind events — defined in that study as occurring if the overnight maximum wind 
(averaged over ten minutes) at Wattle Park was greater than that at Mount Lofty. Grace and Hol- 
ton (1988) showed that the gully-wind often had an easily identified onset and cessation. Namely, 
the wind regime switched sharply from light and variable to fast and fairly steady in direction, or 
vice versa. For Wattle Park, the average onset was at 1815 CSST with the average cessation time 
being 1045 CSST. The mean maximum winds at Wattle Park and Mount Lofty were 10.2 m/s and 
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Fig. 13 Temperature profiles for Adelaide Airport and Mount Gambier at 
0900 CSST on 7 February 1985. ML denotes Mount Lofty summit. 


7.1 m/s, respectively, whereas the mean geostrophic wind was 6.8 m/s. Nineteen out of the thirty 
gully-wind events were associated with moderate to severe turbulence as reported by aircraft op- 
erating out of Adelaide Airport. 


For 29 of the 30 events, low-level winds from Adelaide Airport were used to determine the 
altitude of any low level easterly jet. The wind reports have a nominal vertical resolution of 150 m 
based on an assumed balloon ascent rate of 300 m/min and '/, minute readings. Marsh (1987) 
checked all February 1985 ascents and found that at 0300 CSST the average ascent rate during 
the first few minutes was 310 m/min but with extremes of 240 and 450 m/min. On this basis, Grace 
and Holton (1988) showed that a typical altitude of the low-level jet would have been about 300 
to 400 m. About half of these events had westerly 1000-m winds, and about a quarter had westerly 
or light (less than 3 m/s) winds at 600 m. These results emphasise the shallow nature of the gully- 
wind. 


The author has examined many Adelaide Airport vertical temperature profiles for 2100 
CSST on the night of significant gully-winds. There is no obvious general feature: this is probably 
because on many occasions 2100 CSST is too early and on others, the balloon ascent could be 
either well within the downslope wind layer or in the rebound zone. The nature of this rebound 
zone will be clarified later. Grace and Holton (1988) found that the presence of a subsidence in- 
version on the Adelaide Airport sounding was a poor indicator of a subsequent or simultaneous 
gully-wind. Nevertheless, the most damaging gully-wind events appear to be associated with a 
marked subsidence inversion or a post-frontal inversion. 
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3.2 The case study of 6/7 February 1985 


This subsection is based on the findings of Marsh (1987). On the nights of 6/7 and 7/8 of February 
1985, gully-winds damaged trees and property in many foothill suburbs. Whilst winds were 
strong enough to fell large trees at locations 5 km from the crest, winds at Kent Town, 10 km from 
the crest, were light and variable. Gusts up to 23 m/s occurred at Waite Institute, although about 
40 km further south, at a wind energy monitoring site, gusts to 30 m/s were recorded. At Adelaide 
Airport the maximum wind at 0900 CSST was 110/22 m/s at 400 m AMSL. 


Synoptic charts for 0900 CSST 7 February are shown in Figs. 11 and 12, it should be noted 
that the surface chart therein is very similar to that in Fig. 10. Note that at 0900 CSST the wind at 
1500 m was 100°/02 m/s and at 2100 m was 230°/05 m/s, with westerlies aloft; which again em- 
phasises the shallow nature of the gully-wind. On a synoptic scale, the easterly flow over the 
Mount Lofty ranges ceases at about three times the height of those ranges. 


Vertical temperature profiles for 09900 CSST at Adelaide and Mount Gambier are compared 
at Fig. 13. Immediately noteworthy are the 1400 m subsidence inversion at Mount Gambier and 
the inversion layer between 600 and 1200 m at Adelaide. It is likely that the Mount Gambier trace 
was representative of the stream in general, or possibly a little cooler (being further south), but 
that descent of the sharp inversion occurred over the ranges and to their lee. Shear-induced mixing 
would then act to ‘smear out’ such a sharp inversion. 


3.3 The case study of 25/26 February 1991 


This case is essentially similar to, although not as pronounced as, that in section 3.2. The 0900 
CSST 26 February MSL chart with pertinent 600 m winds is shown at Fig. 14. Contemporaneous 
temperature traces from Adelaide and Mount Gambier Airports in Fig. 15 show a post-trough in- 
version indicative of a cool change penetrating into the southern part of South Australia. The 
air mass which produced the cool change at Mount Gambier is obviously cooler and deeper than 
at Adelaide. 


At 0800 CSST an aircraft at Cleve, about 200 km downstream of Adelaide, reported an in- 
version at ‘1500 to 2000 feet AMSL’, which would be about 550 m AMSL. At Port Stanvac, south 
of Adelaide, the mean wind at 0750 CSST was 14 m/s from the southeast. 


At 0600 CSST, the conditions recorded by the Mount Lofty Automatic Weather Station 
(AWS) were an easterly wind of 8m/s with temperature of 22.5° and dew-point of 2.2°C; yet at 
0728 CSST, they had changed dramatically — to an easterly of 4 m/s with temperature of 12.3° 
and dew-point of 12.0°C! This implies that a 10°C inversion rose through the AWS at 736 m be- 
tween 0600 and 0728 CSST. (This ‘convenient’ occurrence was brought to the author’s attention 
by his colleague, Mr Ian Holton). On a few apparently similar occasions, the Mount Lofty AWS 
records have revealed rapid fluctuations of up to 10° in both dry bulb and dew-point temperatures, 
though at Mount Crawford, at a lower altitude to the north, the sampled air was continuously cool 
and moist. It is evident that on 25/26 February 1991, a reservoir of cold air gradually increasing 
in depth formed a decoupled single-layer flowing across the ranges — very much as proposed 
by Kiittner over 50 years ago. 
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Fig. 14 Charts for 0900 CSST 26 February 1991: (a) gradient level and 
(b) mean sea level pressure chart. 


3.4 The case study of 6/7 February 1972 


On a night without any obvious strong synoptic stream (see Fig. 16), gully-winds were recorded 
by an experimental acoustic sounder at Salisbury by Mahoney et al. (1973). The Salisbury data 
for the four hours from 2200 CSST on 6 February 1972 are presented at Fig. 17 and show the 
height of the turbulent-layer, the vertical velocity at 100 m, and the 100-m layer mean wind 
speed. A measure of consistency is provided by the balloon-derived mean wind to 300 m (the 
lowest available) at 2100 CSST at Adelaide Airport, which was 127°/12 m/s. From the aerolog- 
ical profile at 2100 CSST at Adelaide Airport (Fig. 18), it is seen that the height of the base of 
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Temperature profiles for Adelaide Airport and Mount Gambier at 0900 CSST 
on 26 February 1991. When used in conjunction with Fig. 14, these show that 
there is a shallow layer of cool air crossing the ranges from the east to southeast. 














Fig. 16 Mean sea level pressure chart for 2100 CSST 6 February 1972. 
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Fig. 17. Acoustic-sounding record of turbulent-layer depth, the vertical velocity at 100 m 
and the 100 m-layer mean wind for the period 2200 on 6 February to 0200 CSST 
7 February 1972 at Salisbury, South Australia. 
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Fig. 18 Vertical temperature profile for 2100 CSST 6 February 1972 at Adelaide Airport. 
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Fig.19 Mean sea level pressure charts for 0900 CSST on; (a) 3 February 1988 and 
(b) 4 February 1988. There is weak gradient flow for (a) but a shallow cool- 
air intrusion occurs in the wake of the (dotted) trough in (b). 


the inversion (210 m) was consistent with the Salisbury acoustic sounder data relating to the 
height of the turbulent-layer (about 200 m). The inversion itself was only 10 to 20 m thick. 


Prior to the jump at about 0050 the layer depth was about 200 m. It is less clear what the value 
after the jump was, but a reasonable value would have been about 350 m (Fig. 17). (More details 
of this event are provided by Grace and Holton (1990), of which a copy is included with this the- 
sis.) The remarkable feature is the shallowness of the layer. Corroboration that this is common 
comes from Mr Alan Shore, a hot-air balloon pilot who operates in the Barossa Valley mostly dur- 
ing early morning and late afternoon. (Mr Shore lives near Bethany in the Barossa Valley and 
maintained an anemometer at his property for the author.) According to him, the air aloft during 
gully-winds is calm, yet the surface wind is strong, gusty and noisy. The transition is usually quite 
sharp and occurs at about 200 feet (60 m) above ground level. 
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3.5 Field observations 


3.5.1 Outline of measurements. 


At Adelaide the most common gully-wind direction is from the east-southeast presumably because 
the ranges nearby are oriented north-northeast/south-southwest so that the gully-wind is generally 
perpendicular to the line of the ranges. To utilise Mount Lofty AWS, a 110°—290° transect was se- 
lected as a reference cross-range line. Nairne was selected as a windward aerological and surface 
observation site. Other sites marked in Figs 20-23 were fitted with available instruments. 


On the nights of 2/3 and 3/4 February 1988, designated as nights ‘A’ and ‘B’ respectively, 
Adelaide Airport meteorological radar staff carried out three-hourly sonde and radar-wind flights 
and, at Nairne, Flinders University and Bureau of Meteorology staff carried out three-hourly ra- 
diosonde and visually-tracked balloon flights. 


At approximately coincident times, the Flinders University ‘GROB’ aircraft flew out and 
back along the transect, with several shorter sorties to the lee of Mount Lofty. Because of Air 
Traffic Control requirements, the GROB was unable to fly sufficiently low over the lee foothills, 
and passed about a kilometre south of Mount Lofty (over a well-lit freeway interchange). Instru- 
mentation available on the GROB allowed the acquisition at 13 Hz of air speed, temperature, hu- 
midity, pressure, altitude, and radar height above ground. Navigation and position fixes were 
coordinated by a co-pilot. Horizontal wind was measured, but these measurements were later dis- 
carded as unreliable (the aircraft’s position finding equipment has since been upgraded to the sat- 
ellite-based Global Positioning System (GPS), which obviates this problem). 


Generally, the air/ground coordination was such that the aircraft flew over the balloon sites 
within a few minutes of each balloon release. This high degree of coordination was mainly due to 
the skilled flight-planning of the pilots, for radio transmissions from the aircraft had to be mini- 
mised, as this ‘spiked’ the data from most of the aircraft’s sensors. 


At 2330 CSST on the second night, an additional sonde was released from near Mount Lofty 
AWS (actually about 500 m south of the AWS in order to avoid trees, transmission towers and 
wires), but no wind flight was done. Some other unsuccessful attempts were made. 


Other than the aircraft data, all data were processed using standard techniques and proce- 
dures. The aircra,. data were managed by the pilots, analysed with the package described in Hack- 
er and Schwerdtfeger (1988) and presented as a file containing directly-measured and derived 
parameters with x, y,z,t coordinates. Consistency was found between aircraft and balloon data with 
temperature and dew-point data from various sources differing by less than '/, to 1°C. 


3.5.2 Analysis of study nights A and B. 


Surface charts for 0900 CSST (at the ‘end’ of each night) are shown at Fig. 19. For each three- 
hourly set of data, cross-sections of virtual potential temperature,6,, and specific humidity were 
analysed. These are shown as a series of panels for both nights at Figs 20 to 23. Some subjectivity 
was obviously necessary in the analyses, especially over the lee foothills. The balloon winds were 
also superimposed on these cross-sections. To aid in assimilating the data, the radiosonde flights 
are also presented on conventional aerological diagrams in Fig. 24. These emphasise the strength 
of the inversions on these occasions. 
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Fig.20 Three-hourly isentropic cross-sections (using virtual potential temperatures in 
“C) for night of 2/3 February 1988 along the true bearing 110°-290° through 
Mount Lofty. Faint lines show aircraft tracks and balloon trajectories. Offset 
winds are for Adelaide Airport (left) and Nairne (right). 
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Fig.21 As for Fig. 20 but for specific humidity in gm/kgm. 
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Fig.22 As for Fig. 20 but for the night of 3/4 February 1988. 
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Fig. 23 As for Fig. 22 but for specific humidity in gm/kgm. 
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Fig. 24 
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25 
Temperature 


Three-hourly temperature profiles at Nairne (solid) and Adelaide Airport (broken) 
on the nights of the 2/3 and 3/4 February 1988 (identified as nights A and B, respec- 
tively). At 2330 CSST on the second night a sounding was completed at the crest of 


the ranges. The superadiabatic layers on the B3 profile have been double checked 
and therefore appear to be real. 
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The striking features were that above 1200 m there was apparently little disturbance. After 
either 1800 or 2100 CSST decoupling was generally evident between the stable lower layer and 
the fairly unstable upper layer. The top streamline (usually 8, = 27°C) of the stable layer was typ- 
ically 1000 m above ground to windward. Near the ridge line, it was drawn down rapidly to a min- 
imum of about 300 m over the lee foothills and then passed through a jump-like feature and exited 
at an altitude of 600 m. Wattle Park (10-minute mean) winds varied from 3 to 15 m/s indicating 
that the jump zone was certainly in the vicinity. The specific-humidity analyses and provide a con- 
sistency check with the isentropic analyses. 


Directly in the lee of Mount Lofty itself, e.g. at Waite Institute, Kent Town and at Adelaide 
Airport, the surface winds were mostly only moderate. However, sites at Port Stanvac and Beth- 
any recorded lengthy periods (several hours) of 12-15 m/s and 15 m/s mean winds, respectively. 
At Wattle Park, occasional periods of 15m/s mean wind were recorded. From the author’s friends 
and colleagues, there were many reports that strong gusty gully-winds occurred in the southern 
and northern foothill suburbs. Apparently much of the flow is around the Mount Lofty summit 
and the final panels of Fig. 22 show this effect. This type of effect was noted before in section 3.3. 


3.6 The definition, description and use of AT 


In analyses of actual gully-wind events in the Adelaide region, one of the most important param- 
eters turned out to be the temperature difference, AT, which is defined to be the environmental 
temperature at 736 m above Adelaide Airport minus the temperature at the Mount Lofty AWS 
(which is of course also situated at 736 m AMSL). Marsh (1987) used 0900 CSST radiosonde data 
from the months of December 1985 to April 1986 and determined the average AT, classified ac- 
cording to wind direction at Adelaide Airport (See Fig. 25). It is strikingly evident that in easterly 
streams, the average AT is about 4°C, whereas in others, it is close to zero. 


A somewhat similar relationship is shown by a scatterplot matrix of AT, geostrophic wind 
and peak gust speed at Port Stanvac. The first two variables were determined at 2100 CSST and 
the peak gust used was that occurring within one hour of that time. The scatterplot matrix (Fig. 26) 
shows the dependence of gust speed on AT and on geostrophic wind but little dependence of AT 
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Fig.25 Average AT as a function of Adelaide Airport wind direction at 
736 m, above MSL (after Marsh 1987). 
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Fig.26 Scatterplot matrix of peak gust at Port Stanvac (between 1000 and 1200 
CSST), AT and geostrophic wind speed at 1100 CSST. 


on geostrophic wind. So AT and geostrophic wind are quasi-independent predictors of gust speed. 
It is relatively easy to devise a best-fit quadratic plot of gust speed as a function of AT and geos- 
trophic wind, although this is not shown. 


The greater the amplitude of the forced lee-wave trough then the greater will be AT. In this 
sense, AT may be regarded as a proxy measure (albeit very crude) of the amplitude of the trough, 
and hence of the downslope wind speeds. Grace and Holton (1988) investigated this effect and 
found that, on average, gully-wind speed is a function of AT. Holton and Grace (1989a,b) were 
then able to produce gully-wind nowcasting aids. 


3.7 Observations of a hydraulic jump in a fog bank 


It is well known in hydraulics that transitional flow may occur due to a lateral constriction instead 
of a vertical obstacle (or, indeed, due to both together). An interesting example of this occurred 
in Backstairs Passage at the southern end of the Mount Lofty Ranges on 2 May 1988. More detail 
of this instance is provided in Grace (1991), of which a copy is provided with this thesis. 


In this instance, a shallow cool, moist layer of fog of the order of 200 m deep transi- 
tioned from subcritical to supercritical flow, and eventually back to subcritical via a well de- 
fined hydraulic jump. To the author’s knowledge, this is the only published photograph to 
show clearly an atmospheric hydraulic jump caused by a constriction. One of the photographs 
from Grace (1991) is shown in Fig. 27. Such striking evidence of this type of event leads to 
the conjecture that the actual gully-winds are induced by a combination of lateral constriction 
and vertical obstacle. 
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CHAPTER 4 
A QUASI 11/,-LAYER MODEL (Q17/,) 


4.1 Development 


4.1.1 Assumptions and nomenclature 


Hydraulic models with a single homogeneous layer are naturally termed single-layer models. 
Models with a shallow homogeneous layer overlain by an infinitely deep homogeneous layer are 
formally equivalent to the single-layer models (through the Boussinesq approximation) and so 
are often referred to as single-layer models, even though there are two separate layers in the mod- 
el. To avoid confusion Baines (1987) has used the term ‘one and a half layer’ with the ‘half’ ap- 
plying to the deep passive layer, and this terminology is followed in this paper. 


Consider steady-state 2-dimensional, infinitely-deep, continuously-stratified flow across a 
mountain range as in Fig. 28. The flow is from left to right. The terrain height, h(x), reaches a 
maximum of H at x = 0 and is exactly zero for x 2 O. Far upstream h is effectively zero, and the 
flow velocity, uy, and stability, Ny, are constant with height. The streamline, or potential temper- 
ature surface, 6, originating at level, d, splits on the windward approach to the mountain with 
the lower branch descending rapidly. Above dj), the flow is assumed to be unchanged. Between 
the split streamlines the air is turbulent and well mixed with constant potential temperature, 0), 
but with no mean motion. Other than at the jump, the flow in the layer below the lower split 
streamline is regarded as smooth, non-dissipative, steady, hydrostatic and Boussinesq. 


UNDISTURBED FLOW 


ii 0————— 
8 


0 WELL MIXED 
REGION 














Fig.28 Steady 2D flow over a mountain with upstream static stability Ny, and velocity, uy. 
The streamline initially at altitude d, splits, so that above d, the flow is unchanged. 
Below the lower split streamline the flow is essentially hydraulic, and between the 
split steamlines the air is well-mixed and turbulent but with no mean motion. Sta- 
bility and speed of the layer beneath the lower split streamline (8,) are constant 
with height. 
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The continuously-stratified lower layer is assumed to be approximated by an equivalent ho- 
mogenous layer (of constant density and with velocity independent of height — hence the term 
‘quasi 1'/,-layer’, or Q11/2). No a priori justification for this assumption is attempted other than 
observing that it must be plausible in view of the conclusions of DK and BP; however, it is an- 
ticipated that the realism of the model results will support its validity. Within the hydraulically- 
active, lower layer, the streamlines will remain equidistant and therefore both u and N remain in- 
dependent of height. The infinitely-deep stratified layer above d, is ‘replaced’ by an infinitely 
deep layer of constant potential temperature, 8). It is reasonable to do this since it has already 
been assumed that there are no changes to the flow above d,. A hydraulic jump is assumed to 
occur in the lee of the mountain. Subscripts apply to the active layer in such a way that ‘0’ = far 
upstream, ‘c’ = crest, ‘1’ and ‘2’ = immediately before and after the jump, respectively, and 
“T’= transition threshold. 


4.1.2 Altitude of the splitting streamline 


For an isentropic layer of velocity, u, and of depth, d, overlaid by a deeper isentropic layer, the 
shallow-water gravity-wave speed, c, is given by 


c=utrjgd (11) 


58 
where g’ = g— (12) 
85 


and g is acceleration of gravity, 6, is the upper-layer temperature and 86 is the temperature dif- 
ference between the layers. For a layer of continuously stratified fluid of depth, d, and constant 
Brunt-Vaisala frequency, N, overlaid by an infinitely-deep isentropic layer with temperature 
equal to that at the top of the lower layer, it is known that the gravity waves have a speed given by 


2Nd 


2 yj 
ee eS CA=- Da 


where 7 = 1,2,3,0... (13) 


Equating Eqn 11 and 13 gives an equivalent modified gravity for a continuously stratified layer 
of 


4N°d 
= ——, (14) 
(2n-1) 0 
The layer Froude Number, F, is commonly defined as 
F=s— (15) 





\gid 
although there are several other definitions in the literature. Combining Eqns 14 and 15, we there- 
fore find that 


_ (Qn-1) tu | 
~ 2Nd 


For a shallow homogenous layer overlaid by a deep layer, the condition for transitional flow 
(having F < 1 on the windward side and F > 1 on the leeward) is known to be 


F (16) 
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17 

F 3m? +2(1-7) = 0 ay) 
dy 

(A derivation of Eqn 17 is included at Appendix A.) Parameters d, H and A are non-dimensional- 

ised by setting d=d/ and H= HI where / = N,/uy. (Generally carats will indicate non-dimensional 

values.) The Q1'/2 (Quasi 1'/2-layer) model then has as the condition for transitional flow 


2 * F 
oa (a) +2(1-F) =o. (18) 
2do 2do do 


Equation 18 shows that do is an analytic function of A and nand the graph of this function in Fig. 
29 shows that do increases with H and n. For later discussion the ratio do/H (which is of course 
identical to dp/H) is also plotted thereon. H has only been graphed over the domain (0,1) although 
there is no upper limit imposed by the model. For n=1 (this effectively confines the discussion to 
the lower troposphere), do can be well approximated (using a best-fit power curve) by the explicit 
expression 





i 0.633 1 

dy = 2.278 +5 a) 
while the more general Eqn 18 has a family of solutions for dy with a quasi-period in do of 7 (this 
compares to S85 which has multiple solutions with period 27). A more convenient but less accu- 
rate form of Eqn 19 is 


dy=23H + (20) 


and so 
Uo Ug 
dy =2.3H™ + NW. +m™— (21) 


showing that d, increases with H and u, but decreases with Ny. That is, dy is decreased for a low 
obstacle or if the stream is light and/or stable. A typical value of / is 103m! (with Ny = 10°75"! 
and uy = 10 m/s). So for a 1 km high mountain, dy ~ 4 km and, for a 500 m mountain, d, ~ 3 km. 
The ratio, do/H, or d/H, is about five for H ~ 0.7. 


4.1.3 Critical mountain height 


For continuously stratified flow there is a critical or threshold value (7) such that if A > Ay then 
the flow will adopt the transitional steady-state configuration of Fig. 28. In Chapter 2 it was noted 
that the value of Ar is 0.85 for symmetrical bell-shaped hills (Miles and Huppert 1969; Lilly and 
Klemp 1979). For asymmetric bell-shaped obstacles, Baines and Granek (1990) have shown that 
0.5 < Ar 1: with the lower (higher) values applying as the approach (lee) slope decreases. 
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Fig.29 — Graphical solution of Eqn 18 showing do and higher order solutions as a 
function of H. Also shown is the ratio d/H for n= 1. 


4.1.4 Leeward velocity 


Application of Bernoulli’s equation along the uppermost streamline of the active layer (implicitly 
assuming constant pressure along this streamline) prior to the jump gives 


die ias uo +d (22) 
2g ~ 2g 


and, on substitution for g’ from Eqn 14 (with n = /), this yields 








Tu Tuo 
+dt+th = —.—+d) (23) 
8No do 8No do 
= ik 


where k,, say, is a constant. By continuity it follows that the volume flux (ud) is constant so that 
after some manipulation, Eqn 23 reduces to 


2 


0 
8rd 





tdth = ky. (24) 


Equation 24 is a cubic in d and has two real positive roots which correspond to the sub and 
supercritical flow. Thus, the altitude (d+h) of the uppermost streamline is known at any point 
and, since the layer streamlines maintain their initial relative spacing, the subsequent position of 
a streamline passing through any given height far upstream is determinable. 
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Assuming that no rebound, or jump, occurs before x = , then, setting d = d, and noting 
that h = 0, Eqn 24 yields 


life yv a@ irnry 
ee |" = bala (25) 
2104) do ~\od 


which shows that di, being closely related to do. is also a function of (see Fig. 30). The so- 
called speed-up ratio, §, defined as 





1 (26) 


increases with A (Fig. 31). Figure 31 shows that the leeward velocity is typically 2 to 3 times the 
upstream velocity! Later, it will be shown that $ can be even greater, with favourable topography. 
Under the circumstances, it is not surprising that downslope wind gusts of up to 60 m/s have been 
recorded (Durran, 1986b). Since $ is evidently well approximated by 


ys) 
$=2H +1 (27) 
and therefore 
Ww 
U; = Ut 2uglt (28) 


= Uy t2No Hug’ - 
Thus, for a high-drag regime, u, increases with each of No, u) and H. For example, with uy 
= 20 m/s, Ny = 0.01s"! and h = 1.2 km, then u, = 48 m/s. These are plausible values for the well- 
known chinook windstorms of Boulder, Colorado. 




















Fig.30 Plots of diy dh dy dp» dy. and dis against A. 
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Fig. 31 Plots of jump strength 8, wind-speed amplification 3 and 
turbulence intensity j against H. 


4.1.5 The jump 


If it is assumed that the jump is stationary at x = a and that the leeward and windward plains are 
at the same altitude, then the characteristics of the jump in terms of the upstream parameters are 
easily determined. The basis for the assumption that the jump associated with self-initiating tran- 
sitional flow is located at the point x = is that that point is the first leeward point where h = 0. 
The strongest possible stationary jump occurs at this point: stronger jumps propagate down- 
stream; weaker jumps are stationary on the leeward slopes. Thus this discussion on jumps in 
steady-state flow configuration relates only to the strongest (that is, biggest) jumps. This restric- 
tion is applicable only to the results of subsections 4.1.5 and 4.1.6. From conventional single-layer 
theory it is known that 


= 5A1+8F; -4 G2) 


a) S 
as) 


where F; is the Froude Number immediately before the jump. A derivation of Eqn 29 is provided 
in Appendix B. Adapting Eqn 29 to the Q1'/) model, we then find (after some manipulation) that 


(30) 





Thus, d> is seen to be a function of H and is also graphed in Fig. 30. Note that d> is always less 
than dp, because energy is lost in the jump. Following Lighthill (1978), the strength of the jump, 
B, is defined by 
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d, 7 (31) 


and B is therefore also a function of A (as shown in Fig. 31) and can be approximated by 


b= 1sJff (32) 


showing that it increases with either of N, or H but decreases with uy. Of course, if the jump is 
positioned at some point up the hill, the energy losses, $ and B are decreased. The energy is dis- 
sipated by either dispersive short waves downstream from the jump or by turbulence at the jump 
face. It is of course the turbulent type of jump which is of especial relevance to aircraft. For 
B > 0.75 turbulence is the main agency of energy loss for water flow. The corresponding value 
of B for two-layer flow is probably significantly greater but not yet established accurately (Rott- 
man and Simpson 1989). It is speculated here that for stratified flow the threshold value of 

B is 1.38 which from Eqn 32 would imply that H = 0.85. That is, possibly A turns out to be 
0.85 because B, is 1.38. Once Bis = 1.38 then sufficient turbulent mixing would sustain the de- 
coupling between the hydraulically-active lower layer and the well-mixed passive layer (passive 
in the hydraulic sense). 


The energy dissipation rate per unit breadth of the jump, e, at either type of jump is given by 
Lighthill (1978) as 


aady 


e@ = pug 4d 
2 


(33) 


where is density. For a turbulent jump, the dissipation rate per square metre of jump face is de- 
fined here as the turbulent intensity, 7 such that 


j=. (34) 
Therefore, after some algebraic manipulation, we find that 


pug do (dy -d,)* 
j= 





5 (35) 
nd,d, 
Define normalised turbulent intensity,7, by 
je (36) 
Puo 
then 
Poe 
> _ d)(d2-d)) 
2 G7) 
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j is a function of H (see Fig. 31) and is approximated by 


j= 148 , (38) 
and thus, 
j = 1.4 (NoHuy)*”. (39) 


It is therefore expected that as any of No, H, or uy increase, then the intensity of turbulence at 
the jump face increases; this result is certainly intuitively reasonable. Of course, if the jump re- 
mains undular then it is the amplitude of the undularities which will increase. 


4.1.6 Wave drag 


It is convenient to define a normalised wave drag, 7, as 


rNo 


Ml 


ra 


; (40) 
p uo 


where r denotes wave drag. Other common terms used to describe this quantity are mountain drag, 
mountain wave drag, pressure drag and gravity-wave drag. The wave drag, r,, associated with lin- 
ear waves over a symmetric mountain is given by (e.g. Gill 1982) 


| 2 
r= zPmNowoll (41) 


so that 


(42) 


Assuming that any jump on the lee side is stationary and located at x > o, then, for the Q1!/, mod- 
el, the normalised drag, 7, is (see Appendix C) 





2do ] A 
p= 7 (1-3 |-ae@-0 (43) 


iat 


The graphical comparison in Fig. 32 shows that this drag is more than three times the linear drag. 
Values from BP and DK are also plotted, and it is clear that these compare quite well with the 


1.95 
Q1!/, model drag. Since * can be approximated closely by the curve 26/7 and the latter expo- 
nent can be rounded to 2.0, it follows that 

r=3.3r,. (44) 


Thus 7 increases with any of Ny, uy or H. 
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Fig.32 | Nondimensional drag, 7, for S85 (labelled S), Q1!/, model (labelled Q) and linear waves 
compared as functions of H with spot values from numerical experiments by Durran 
and K lemp (1987) (triangles) and Bacmeister and Pierrehumbert (1988) (circles). 








7 wl 

Fig.33. Model solutions of do for Q1!/, and $85/SS87 as a function of H with spot val- 
ues from tank experiments (unfilled circles) and numerical experiments (filled). The 
spot values and S85/SS87 curve are from Rottman and Smith (1989). 
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4.2 Comparison with Smith’s model 


The main difference between the assumptions of S85 and Q1!/, is that S85 does not assume that 
the velocity in the hydraulically active layer is constant with height. Partly because of this, S85 
makes no comment regarding any eventual rebound of the accelerated flow. The predictions of 
S85 and Q1!/, for both do and dyare fairly close (see Fig. 30) with the Q1'/, values being up to 
10-20% lower. From data provided in S85, it is easily shown that, at least to good approximation, 


dys =nH +- (45) 


where the subscript ‘s’ denotes S85 model. This is similar to the Q1!/, approximation at Eqn 20 


and its simple form suggests that an alternative derivation to that given by S85 and SS87 is pos- 
sible. However, this has not been attempted herein. 


Values for d, from the numerical models are known to accord fairly well with S85 values, 
which in turn are about 10-15% higher than the Q1!/, values (see Fig. 30). However, the vertical 
resolution in the models is usually about 200 to 350 m and this is too coarse to determine clearly 
whether S85 or Q1!/5 gives the better estimate. Rottman and Smith (1989) compared predictions 
of dy from the model of $85 and SS87 with their tank experiment results and with some numer- 
ical model results from the literature. Figure 33 shows the S85 do value, the tank results, the nu- 
merical model results (adapted from Rottman and Smith) and the Q1!/, model predictions for do. 
Obviously, the Q1!/, values compare extremely well with the experimental data; but S85 ones 
are equally as good. 


_ Interestingly, the speed-up factor, 3, is effectively identical for both models. For example, 
at H = 1 the S85 value is 3.017 (from SS87) compared to 3.022 for the Q1!/, model. This im- 
plies that the layer mean wind $uy over the lee foot of the hill is virtually identical in each case. 


S85 has multiple higher-order solutions of the form do +2nn, where n=0,1,2,..., and these 
appear to correspond approximately to the well-elevated but narrow bands of environmental crit- 
ical-layer positions associated with a weak tendency to initiate transitional flow. Q1!/, has mul- 
tiple higher-order solutions of the very approximate form dy +m where n = 0,1,2,... Compared 
to S85, Q1!/, has too many solutions in as much as those for odd n are invalid. More importantly, 
this points to a specific area in which to test Q1'/, against S85 using numerical and tank models. 
Even so, it seems that at least for the Adelaide gully-winds where, as shown in Chapter 3, dj is 
less than about 2h, then the higher terms are not of practical significance. 


For S85, the drag, rs, is given by 


r= oN (dy-d,)* (46) 


or, equivalently, 


>. & ‘Fh 
t, = 2 (dp=a) , = (do- 5) (47) 
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where the subscript S denotes S85 values. 7, is graphed as a function of A in Fig. 32, where it is 
seen to be significantly greater than that of the Q1!/, model. In comparison with the numerical 
model results, S85 overpredicts the drag by a factor of 2, and this indicates, through Eqn 47, that 
it overestimates dy by about 15%. As mentioned, the vertical resolution of the numerical models 
is such that this is not easily verifiable, at least from results in the literature. 





Fig.34 The defining sketch for M. The topography has been modified in such a way that the 
lee plains are lower than the windward ones. 


4.3 Topographic effects 


Now let us consider the topographic case in which the lee plains are lower than the windward 
plains by an amount MH, where M is defined as ( h_,,—h, ,,) /H, as in Fig. 34. In such a case, 
h=-MH and d= d, at x = Q, so that Eqn 23 yields 


2. 
~ 
sae 


8No dy 





+d,—-MH = ky (48) 


so in similar fashion to the derivation of Eqn 25, we can find that 


17m) ad MA i7ny 
i = Gales (49) 
2d;) dy a oa 


In this case, d, and hence, 5, B and 7 are functions of (M, A) (see Fig. 35). In particular it is ev- 
ident that w,, B , and j increase with M, that is, with a ‘dropped’ lee plain. This result is intuitively 
reasonable and consistent with the occurrence of chinook windstorms along the Front Range in 
Colorado. 
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Fig. 35 
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Plots of (a) jump strength, B, (b) wind-speed amplification $ and (c) turbulence 
intensity 7 against H for various values of M. 
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Now consider a symmetric range comprising a string of linked conical hills as in Fig. 36. 
The relevant feature of cross-range flow is that there are two extreme possibilities, namely, val- 
ley-valley (VV) flow or spur-spur (SS) flow. For valley-valley flow, Arakawa (1969) extended 
the conventional hydraulic theory and this extension may be applied to the Q1!/, model. Conse- 
quently, Eqn 18 (with n = 1) becomes 


2 23 mM 
ea =o +2(1-2) = (50) 
2do 2d do 


where y= m_./m, and is greater than one, and m- is the (two dimensional) mass flux at the crest. 
The parameter, Y, may be regarded as a valley funnelling factor. Sufficiently far downstream, the 
flux, m, again equals m,. Arakawa has shown that the windward convergence/leeward divergence 
in such flow implies that the valleys experience more violent downslope winds than the spurs. 
More recently. Pettre (1982) has used the concept to model successfully the Rhone Valley winds, 
and Grace (1991) used a similar idea to model transitional flow through Backstairs Passage (see 
section 3.7). Strictly speaking, Arakawa’s model applies to an homogeneous layer bounded by 
vertical valley walls. However, it seems reasonable to apply it qualitatively to the Q1'/, model. 
The concept may be extended to SS flow (y < 1) where windward divergence/leeward conver- 
gence is expected. It is reasonable to assume that all other flows normal to the range will lie be- 
tween the VV and SS cases. Equation 50 shows dp to be a function of Hand y. Interestingly 
enough, transitional flow may occur with H = 0. The fact that dp increases with y, 1s a result 
consistent with Arakawa’s findings relating to F, in the homogeneous layer case. 














Fig. 36 Plan view of a string of conical hills ‘melded’ together. Heavy solid lines are con- 
tours and dashed lines represent streamlines. Flow along VV represents flow 
through a valley while SS represents flow along the highest path or spur line. 
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If the valley mouth, the jump, and the point where m again equals mz, are all collocated at 
x =a then, as with the lower lee-plain effect, each of ais’ B and 7 may be determined as func- 
tions of FAA and y (see Fig. 37 - but note that d, not shown). Of special interest is the fact that w,, 
B and 7 increase with , that is, with valley-valley flow. This is a similar effect to that found with 
lower leeside plains — again, an intuitively expected result. The occurrence of Antarctic ‘katabat- 
ics’ is also consistent with this result. 


4.4 Critical-layer effects 


4.4.1 Development 


The numerical experiments of DK and BP show that even for subcritical mountains (those with 
Hf < Hy), the high-drag regime may occur if the model atmosphere has an environmental critical- 
layer. For conventional single-layer hydraulic flow, Houghton and Kasahara (1968) (hereafter 
HK) have shown that, in reference to their diagram (see Fig. 38), transitional flow occurs when 
Domain IIa or IIb applies, but does not occur for Domain I. In Domain Ib, propagating jumps 
are expected whereas, in IIa, the jumps are stationary above a point on the lee slope at which h/ 
H corresponds to the dashed isopleth value. An advantage of the Q1!/, model is that this diagram 
may be adapted by assuming that the environmental critical-layer is analogous to the free surface 
in the single-layer model of HK and that the appropriate g’ for the layer is given by Eqn 14 with 
n= 1. By way of precedent, it is noted that Klemp and Durran (1987) successfully applied HK’s 
diagram to the case of a strong inversion layer embedded in an otherwise continuously stratified 
atmosphere in order to demonstrate the hydraulic-like character of the flow simulated by their 
numerical model. They did, however, use a modified gravity calculated by setting the mid-layer 
6 value to the average 8 value for the layer. This led to 


, ee) 
& = 5Nodg (51) 


in contrast to that of the Q1!/, model, in which (from Eqn 14) 


g ~ ined, = = 0.4Nid- (52) 


T 


where d, is the environmental critical-layer height far upstream. The vertical axis of Fig. 38 in 
the HK model is the Froude number of the single-layer far upstream, and this is equivalent to 
"Vonldg in the Q1!/, model. The horizontal axis is the ratio of the mountain height to the single- 
layer depth far upstream and this becomes H/d, or H/dg. The special case dg = dp coincides 
with the jump/no jump boundary between Domains I and IIa. Defining d, to be that value of d; 
for which the jump is stationary at x = & (where h = 0) corresponding to the boundary of Ila/IIb, 
it can be mapped from HK’s diagram (Fig. 38) and may therefore be regarded as a function of 
(Figs 30 and 39). 


4.4.2 Results 


There are three situations for a subcritical mountain height. The first is that there is no transitional 
flow if dg > do. Otherwise, if dg < do then upstream adjustment occurs with a bore propagating 
indefinitely upstream. The values of d;, No, Uy (but notably not g’, which is unchanged) are thus 
effectively modified by the flow itself to set up transitional flow with a leeside hydraulic jump. 
The second situation is that if dp < dg < do, then stationary jumps occur on the lee slope. Thirdly, 
if dg < dp then jumps propagate downstream over the lee plain. This compares well with 
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Fig. 37 Plots of (a) jump strength, B, (b) wind-speed amplification § and (c) turbulence 
intensity j against A for various values of y. 
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Fig. 38 Adaptation of the diagram of Houghton and Kasahara (1968) for transitional and 
nontransitional flow regimes to the Q1!/, model. 











Fig. 39 Comparison of the results of Durran and Klemp (1987) for high drag (filled cir- 
cles) and low drag (hollow circles) events with Q1'/, predictions. 
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the numerical model results. Unfortunately, there is some ambiguity in the precise definition of 
d,;, as we do not know whether it should be taken to be at the top, bottom, or some mid-point of 
an environmental critical-layer. If the bottom is selected, however, there is then reasonable agree- 
ment between the numerical and Q1!/, models. 


4.4.3 Comparison with Durran & Klemps’s numerical results 


The agreement between Q1!/, and other models is illustrated by the plot of high-and low-drag re- 
gimes from DK’s simulations given in Fig. 39. There is a one-to-one correspondence between 
high-drag regimes and downslope wind events (DK). The Q1!/, model predicts that if A<4Hr 
then the curve dg = do separates the two drag regimes (low drag above, high drag below), and 
assumes that if H > Hr then high drag occurs regardless of d,;. From this point of view, the 
Q1'/, predictions are therefore quite good. 


The lowest plot at H = 0.1 showing a low-drag regime is an apparent discrepancy, but is 
arguably a high-drag regime on the grounds that BP’s simulated outcome for H = 0.08 was one 
of high drag. BP suggested that their different result was because their simulated run time was 
significantly greater. In BP’s simulations, an hyperbolic-tangent wind profile was used and so any 
selection of a bottom point in the environmental critical-layer would have been quite arbitrary. A 
detailed comparison with their results has therefore not been attempted herein. The numerical 
models indicate that there are also well-elevated but narrow bands for d; (above 13 km) which 
have a weak tendency to initiate transitional flow. 


4.5 Summary 


An important assumption in the model is that the hydraulically active layer can be approximated 
by a single homogeneous layer; in effect, the flow velocity is assumed uniform with height. This 
assumption allows the application of a comprehensive body of well-established single-homoge- 
neous-layer hydraulics. It appears justified a posteriori, since the model results compare reason- 
ably well with S85 and quite well with sophisticated numerical models. 


It was shown that as each of N,, H and uy increases, so too does u,. For turbulent jumps, in- 
creasing N,, H or u, increases j, the intensity of turbulence at the jump face. For undular jumps, 
increasing Ny, H, or u, increases the amplitude of the dispersive short wave undularities associated 
with a jump. In a speculative vein, since jump strength, B, was shown to be a monotonically in- 
creasing function of normalised mountain height, A, it appears that increasing A increases the 
likelihood that B will exceed a threshold value, B,, and therefore that the jump will become tur- 
bulent. In this respect, B, could be linked to Hin such a way that self-sustaining decoupling oc- 
curs — a concept similar to that of PC. 


Consideration of the relative elevation of the lee-side plain showed that turbulence intensity, 
the jump strength and uw, increased with a lower lee-side plain. A similar effect was found for fun- 
nelling due to valley-valley flow. It is therefore expected that valley-valley topography with lower 
lee-side plains (relative to the windward plains) will be especially prone to strong downslope 
winds. 


In comparison with $85, Q1!/, underpredicts dp by up to 10-20% although, within resolution 
limitations, both compare well with numerical and tank-model results. Drag estimates from the 
Q1'/, model are 50-60% of the S85 estimates. The Q1!/, drag compares very well with the nu- 
merical models. Compared to $85, Q1!/, has ‘twice too many’ higher-order solutions for dy. In- 
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terestingly, the speed up factor, § = d)/d, is effectively identical for both models, implying the 
same layer-mean speed over the lee plains. 


In the model, a high-drag regime is assumed to self-initiate if H > H;. However, if H < Ar, 
then the model predicts that a high-drag regime is initiated if an environmental critical-layer is 
located below or at d, (that is, dp < dy). Overall, the very simple Q1!/, model therefore yields re- 
sults similar to those of the sophisticated numerical models when simple upstream profiles are 
specified. 


CHAPTER 5 


A MODEL OF GUST PERIODICITY 


5.1 Background 


A feature of downslope winds is their directional steadiness and the gusty nature of the speed. 
A striking exception is the Antarctic ‘katabatic’, which is renowned for its steadiness of both di- 
rection and speed. Boras of Yugoslavia tend to have a quasi-periodicity of 3 to 4 minutes, al- 
though periods of 3 to 11 minutes have been reported. Lidar observations by Nieman et al. 
(1988) of a Boulder windstorm showed a tendency for gusts and surges to occur with periods 
of 3 to 4 minutes and of 14 minutes. Scinocca and Peltier (1989) performed experiments on Pel- 
tier and Clark’s numerical model using a much larger domain than in previous experiments. 
They were able to demonstrate quasi-periodic transience in the wind speed in the 5 to 15 
minute range, and attributed it to the breaking of the topographically forced internal waves. 


5.2 Development 


For stratified fluid the vertically propagating internal gravity waves have a frequency of oscilla- 
tion (in the vertical) of /N/n. Horizontally propagating waves in otherwise still fluid have a 
speed, c, given by 





d 

> = 2N- 53 

c . (53) 

and a frequency of oscillation f (in hertz) of oscillation in the vertical given by 
iB 

ad 3 ai 
_ 2Nd 
oT 


where A is the wavelength. We assume that topographic irregularities induce internal gravity- 
waves in the continuously-stratified fluid. If we further assume that an (unspecified) mechanism 
due to resonance causes the vertically-propagating waves to feed energy to the horizontally- 
propagating waves, then, of all the possible combinations of frequencies, the one where the two 
types of waves have the same frequency is the most important. So suppose that the horizontal 
waves somehow ‘tune’ to the same frequency as the internal waves: after all, they are both os- 
cillating in the vertical; then 


r= Ad. (55) 


Superimposing a flow speed of u, then the local frequency becomes 
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The dual values of f are positive when 
Tu >2Nd 
TU 


anda? 1 (57) 


that is, when F' > | (see Eqn 16); which (as expected) states that waves in the supercritical region 


of the fluid cannot travel upstream. Denoting leeside parameters by the subscript 1, a special case 
of Eqn 56 is 


uy x N, Sy 
~ Gd," In o 
A property of the Q1!'/, model is that 
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5.3 Results and discussion 


From Fig. 31, values of § are in the range 2 to 3 and dj is in the range 2.5 to 4. Choosing as 
typical, § = 2 and dy = 3 and Ny = 0.01s"' then f = (0.56+0.10) 10°*s-! corresponding 
to periods of 3 min and 17 min. 


By using Eqns 20 and 27, namely, 
‘ Tt 
djp=2.3H += (62) 


and 


$= 2H 41 (63) 


itis evident that f is a function of (No, A) which may be expressed as a plot. More conveniently, 
if we define the period by P = f ~! and then plot P against H/ug for selected values of No; we obtain 
the diagram shown in Fig. 40. It is apparent that for H = 0.5 and No= 0.01s"!, the leeside peri- 
odicities expected are 3 and 17 minutes. In view of the periods actually observed, these are very 
reasonable values. It is also apparent that the shorter (longer) period decreases (increases) with 
increasing mountain height. 


Such reasonable results indicate the plausibility, at least, of the hypothesised resonance in- 
teraction, whereby internal gravity-waves, induced by topographic irregularities, feed energy to 
the horizontal waves. The fact that the Antarctic downslope winds are exceptional in their stead- 
iness could be due to the presence of a relatively smooth and friction-free ice surface. 





Period, P {minutes} 

















Fig. 40 Plot of the dual periods in minutes against H/ug and No. 


CHAPTER 6 


THE TIME-DEPENDENT 1'/,-LAYER MODEL (QT) 


6.1 Background 


For continuously-stratified flow with an environmental critical-layer positioned below do, the 
flow is hydraulic in nature (BP). For many downslope winds (e.g. the bora and the Adelaide 
gully-wind especially), the height of the environmental critical-layer, or discontinuity, is well 
below do, so that it is unlikely that PC’s mechanism is important therein. To clarify this point, it 
is noted that Peltier and Clark’s mechanism predicts a decoupling height, do, close to that of 
S85 and Q1'/2, so the predicted value of do/H is about 5 (see Fig. 29). However, observations 
of several events (see Figs 4, 20 and 22) show that the decoupling heights are much lower. 


In these circumstances, it is natural to ask whether the steady-state Q1'/, model can be ex- 
tended to include frictional and Coriolis forces and temporal variations. The answer is in the af- 
firmative and the model developed herein (for convenience hereafter referred to as the QT 
model) is a variant of those presented by Houghton and Kasahara (1968), Parrett and Cullen 
(1984) and Baines and Leonard (1989). In this model, the value of g’ is determined as appropri- 
ate from either Eqn 12 (for a homogeneous layer) or Eqn 14 (for a continuously-stratified lay- 
er). At relatively low resolution, say 100 to 200 grid-points, it is capable of running in an 
interactive mode on a desktop computer. A floppy-disk version suitable for a ‘386’ PC with a 
numeric coprocessor is included with this thesis. 


Typically all three models above consider an obstacle of height H perturbing the otherwise 
uniform flow of a homogeneous layer of initial height dp. The fluid is impulsively started with 
an initially free-standing surface and in all cases (to the author’s knowledge, at least) dg > H. 


The model developed here will be used to investigate two ‘archetypal’ situations: first, a 
geostrophic flow perturbed by topography (Type I) and second, a very shallow layer (where dp 
< H) initially flowing parallel to a mountain range and not in geostrophic balance (Type II). 
The simulations are relevant to the case studies of Chapter 3. 


6.2 Development 


The appropriate starting equations are the shallow-water equations as presented by, for exam- 
ple, Pedlosky (1987) with additional quadratic drag terms. These are the two horizontal momen- 
tum equations and the mass conservation equation, namely: 








CpUu 
u,+uu,+vu,+ gh, + gd,—fv+ q -P,-EV(u,—u) = 0 (64) 
CpUv 
v,+ uv, + vv, + gh, + gd, +fut ae -EV(V,-v) =0 (65) 
d,+ (ud), + (vd), = 0 (66) 
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where (u,v) is the fluid velocity in the lower layer, U is the magnitude of that velocity, (u4,Vv,4) 
is the fluid velocity in the upper layer, V is the magnitude of the shear between the layers, 
h(x,y) is topographic height, d(x,y) is fluid depth, f is the Coriolis force (negative in the South- 
ern Hemisphere), Cp and E are surface and interfacial drag coefficients, and P, and P,, repre- 
sent external pressure-gradient forces in the x and y directions. Subscripts x, y, t denote 
differentiation. The interfacial drag is the drag between the lower and upper layers. It is analo- 
gous to surface drag but unlike surface drag does not necessarily have a retarding effect. In the 
case of a fast-flowing upper layer, for instance, the interfacial drag could act to accelerate the 
lower layer. 


It is these equations, without the drag terms, which are the basis of the well-known con- 
servation of potential vorticity rule and this in turn provides a qualitative check upon the re- 
sults. Before proceeding, the assumptions are made explicit. 


The underlying assumptions are that a layer of fluid (usually thought of as water) is in- 
compressible and of constant density; the layer has an aspect ratio (defined as that of the verti- 
cal to the horizontal scale) very much less than unity. The hydrostatic assumption (that the 
distribution of pressure with respect to height is unaffected by vertical motion in the fluid) is 
also invoked. Usually g is replaced by Eqn 12 (g’ = g60/8,) for the case of a shallow homo- 
geneous layer but, as shown earlier, it would seem reasonable to replace g according to Eqn 14 
(ig = 4N7d/n) in the case of a shallow continuously stratified layer. 


Now consider a mountain range invariant in the north-south direction associated with a 
flow for which uw, v,.and d,, are zero. In this case, the three equations reduce to 








h+d Cott op _ py 0 

+uu.+ +d) —fv+ -Pp- -u) = 

u,tuu,+g(h die =I 7 , (uy — u) (67) 
CpUv 

v,tuv, t+ fut 7d —-P.-EV(v,-u) = 0 (68) 

d,+ (ud), = 0. (69) 


It is always possible to derive an energy equation from the two horizontal momentum equa- 
tions. Since hydraulic jumps with energy losses are anticipated, it is clear that energy will not 
be conserved. For this reason, then, we shall work with the flux form of the momentum equa- 
tions. So we add (Eqn 68) d to (Eqn 69) v, which gives 


(vd), + (uvd), + fud+ CpUv — Fid-EV (wv, -v) =0 (70) 
and we add (Eqn 67) d to (Eqn 69) u, which yields 


(ud), + (wd), + 1. (da), + gdh, —fvd (71) 


+CpyUu-P.d—-EV(u,-—u) =0. 
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The three underpinning equations for the model are now therefore Eqns 69,70 and 71. 
This system of prognostic equations may be regarded as initial-value flux-conservation prob- 
lems and solved numerically to second-order accuracy by using centred time and space differ- 
encing (Press et al. 1986). Parrett and Cullen used this type of technique; specifically, they used 
a variation of a staggered leapfrog scheme. However, their scheme of equations did not incorpo- 
rate drag or external forces: this author found that with values of drag (greater than about 
0.003), their technique was unstable, apparently due to mesh drifting. Baines and Leonard used 
a more sophisticated upwind differencing scheme. The method presented here is quite simple 
and is a variation of the Lax-Wendroff two-step scheme which is described in Press et al. 
(1986). It was found to be stable for the range of parameters tested. A pragmatic feature of this 
model is that at low resolution it is stable and will run in a semi-interactive mode on a desktop 
computer. The author and colleagues (professional weather forecasters) have trialled the soft- 
ware in a subjective fashion. Although ideally it needs accurate upstream data, the input of sub- 
jectively estimated upstream conditions gives results consistent with experience. 


If we define m,n,o,u and A as 
m=ud u=mu (72) 
n=vd d=nu 


ot 


hee , 


then provisional ‘halfway’ values of m, n and d are calculated from: 
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2d 2 = dj, td;—A (ij, ;—m) 


pS 


(75) 
where m., ni, 0, Hi; represent discrete values on a time-space grid. 
Using these provisional values, updated values are calculated from the properly centred expres- 


sions: 
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where the final term in each is a numerical diffusion term. « is the numerical diffusion coeffi- 
cient and is chosen so that 


he < i. (79) 


Chapter 6 D7. 


Finally, it is known that a further stability criterion is 


r (80) 


1 
<. 
(Bias ote V | 
where subscript ‘max’ refers to the maximum value. 
Flow which is initially in synoptic balance (the external pressure-gradient force balanc- 
ing the combined drag and Coriolis force prior to the introduction of topography) has P, and 


P,, determined as 


Cp Upuy — EVy (Ug — Uo) 





P. = —fvgt (81) 
; 0 dy 

Cy Uv, — EV, (V4 — Vo) 
P. = —fuy+ pYo0"o 0\%a— Yo 





dy 


where subscript 0 refers to initial upstream conditions. For flow which is not initially in such a 
balance (e.g. a sea-breeze) then P, and FP, are balanced against the drag in such a way that 


Cp Upttp~ EV, (uy — U9) 








= (82) 
x dy 
CpyUgvy — EVo (V4 — Vo) 
P= d ; 
0 


6.3 Comparison with other models 


Four comparisons were made with results from the model of Houghton and Kasahara (1968), 
without rotation and or friction. With «A/dx = 0.07, the results were practically indistinguish- 
able. 


Parrett and Cullen (1984) used the same topography and initial conditions as Houghton 
and Kasahara, but included unbalanced rotation. Again, the results of their model when com- 
pared with QT model results were practically indistinguishable. A comparison against Baines 
and Leonard’s (1989) two sets of plotted results also yielded very close agreement. These au- 
thors used asymmetric topography and an external force such that the unperturbed flow was in 
geostrophic balance. There is sufficient evidence therefore to accept the numerical veracity of 
the QT model. 


6.4 Simulations using the QT model 


An east-west cross-section (looking south) through the Mount Lofty Ranges was simulated by 
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a H for x >0 (83) 
2 
Li 
a 
h= 2 forx <0 
xe 
l+— 
b 


with H = 500 m, a = 3 km and b = 15 km. Thus the flow is from left to right over an asymmetric 
ramp-shaped hill. Note that easterly and northerly wind components will be positive in this ref- 
erence frame. 


From Stull (1989), typical values of surface drag coefficient in neutral conditions over 
plains are about 0.002: over the sea, Cp is 0.0007. Raupach et al. (1991) have measured drag 
and associated parameters over Australian farmland, forest and city surfaces. According to 
Raupach (personal communication), a plausible value of Cp over the well-treed parts of the 
Mount Lofty Ranges and the built-up areas of the lee plains would be about 0.010. Some of the 
values cited by Stull for coniferous forests are even an order of magnitude greater! At night, de- 
coupling would minimise Cp whereas surface heating during the day would maximise it. 


The left and right end points were chosen sufficiently distant to avoid interference from re- 
flected perturbations, nevertheless, zero-gradient boundary conditions were applied. The ratio 
KA/dx was 0.06, where 5x was between 500 and 700 m. In the results, only a windowed, or 
zoomed, section is shown. 


6.4.1 Balanced synoptic flow — (Type I) 


Experiments with a ‘Type I’ stream in geostrophic balance introduced to topography (with and 
without friction) were conducted. Parameters were chosen to simulate a homogeneous layer in 
synoptic balance beneath a 6°C subsidence inversion at 850 Pa. As ‘reference’ or standard syn- 
optic conditions, uy = 10 m/s, vg = 0, g’ = 0.2 ms, dy = 1500 m, f= -8.3x10°s"!, Cp = 0.002 and 
E=0 were taken. P, and P,, are determined from Eqn 81. Note that, since two homogeneous lay- 
ers were used, this model was essentially a conventional 1'/;-layer model with drag and inertial 
effects included, however. 


The results of a preliminary run (denoted R41) with fand Cp set to zero are shown in Fig. 41. 
(Run numbers and figure numbers correspond). The four panels show the evolution of the flow in 
terms of u,v and d at three-hourly intervals. The features of note are the two hydraulic jumps, the 
maximum wind of about 23 m/s which occurs just upstream of the leeward hydraulic jump, and its 
steady leeward movement. Although not shown, one-hourly ‘snaps’ illustrate this point quite clear- 
ly. Meanwhile, the windward flow speed decreases significantly downstream of the upstream jump 
as it propagates upstream. 


A second run, R42, this time with f = -8.3x10%s"! for latitude 35°South, was performed. 
Compared with R41, the leeward jump accelerates leeward and appears stronger, whereas the 
windward jump decays and decelerates, eventually reaching a steady state (although on these pan- 
els only the region 100 km to windward is shown). The decelerated flow veers and the accelerated 
flow backs — in accord with potential vorticity conservation in the Southern Hemisphere. 
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Figures 41 to 49 — evolution of Type I flow with ‘Witch of Agnesi’ topography given by 
H = 500 m, a = 3 km, b = 15 km, and preconditions of up = 10 m/s, vo = 0 m/s, dy = 1500 m, 
g = 0.2 m/s’, f= -83 x 10°s"!, Cp = 0.002 and E = 0.0 unless noted otherwise. 
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Fig. 41 The case in which f= 0 s"!, Cp = 0, in which transitional flow 
and downslope winds eventuate. 


A third run, R43, which was used as the reference for all subsequent comparisons involv- 
ing synoptically balanced flow, had the same preconditions as R42 except that Cp = 0.002. In 
comparison to R42, the motion of the jumps has been slowed and the maximum wind slightly 
decreased. 


R44, with ug reduced to 5 m/s, exhibits a nontransitional flow. The stronger wind zone is 
confined to the crest region and there is no jump. R45, with dp raised to 2000 m and ug reset to 
10 m/s yields similar results. 


To investigate the effect of inversion strength, g’ was set to 0.1 and 0.3 m/s? in runs R46 
and R46, respectively. R46 has weaker maximum winds, but the jump has propagated further 
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Altitude !km) 








-100 -80 


Fig. 43 








-100 -80 -60 


Distance from crest (kin) 


° faerie eee 


-40 -20 


The ‘reference’ preconditions for Type I flow. 








(m/s) 


wind speeds 


Component 


42 The case in which Cp = 0 but f= -83 x 10s"! In this, the leeward jump 


Wind speed (m/s) 


(km) 











Chapter 6 





(m/s) 


61 





Altitude 








4-—- ne = 1 0 Ajeet 80 
t= 6 br 20 j t-Olr 20 
3 oe eee 10 3 
Pes eae ea tama creel 
4 2 














Qt errr 


-100 -80 -60 -40 





-20 0 -100 -80 -60 


Distance from crest (km) 


Component wind speeds 


(km} 


Altitude 


Fig.44 Non-transitional flow for the reference conditions except that 


u=S5nvs. 











-40 0 




















-80 40 -80 -40 0 40 
Distance from crest (km) 
Fig.45 Non-transitional flow for the reference conditions except that 


dg = 2000 m. 


(mis) 


wind speeds 


Component 


62 Meteorological Study 4] 


leeward. R47 is nontransitional. It is therefore evident that increasing g’ increases the maxi- 
mum wind in the supercritical zone but decreases the likelihood of transitional flow. This latter 
run appears to simulate well the first case event described in section 3.2, where foothills sub- 
urbs experienced damaging winds and yet a few kilometres downstream the winds were light. 
Note, however, that modifying the conditions of R47 by lowering dp to 1000 m causes very 
strong downslope winds to propagate further leeward (though this run is not shown). 


From the above experiments, it can be inferred that increasing up/dp acts to increase the 
likelihood of transitional flow, to increase the maximum wind and to push the leeward jump (if 
any) further leeward, thereby exposing a greater area of the lee plains to downslope winds. 
However, increasing g’ increases the subsequent leeside wind speed maximum but decreases 
the leeward extent of strong winds. 


Numerical modelling by Richard et al. (1989) and by Pitts and Lyons (1990) has shown 
that the inclusion of surface drag acts to delay and to weaken downslope winds. By comparing 
R42 and R43, it can be seen that the effect of surface drag is to weaken the maximum wind at- 
tained and to limit its leeward extent. Not only are these results intuitively reasonable, but they 
are also consistent with those of the sophisticated models. In run R48, the flow is initialised as in 
R43 but the value of Cp in the vicinity of the mountain is greatly increased. Ten kilometres ei- 
ther side of the crest, Cp is set to 0.01, in keeping with Raupach’s estimate, and is then reduced 
to 0.002 with a single cosine step over a further ten kilometres. That is, the perturbation is now 
due to a hill which is rougher than the plains. The dramatic effect of having a relatively rough 
hill is evidenced by the weaker jump, the weaker wind maximum and the lesser leeward extent 
of the jump. This effect was further investigated by means of another run, R49, which was initial- 
ly identical to R48, except that Cp was halved over the entire domain. Here it is seen that the 
wind maximum is greater and that the jump is located further leeward. Thus, surface drag acts to 
inhibit downslope wind maxima and to draw the strong wind zones back toward the crest. 


It was conjectured that friction might account for the known ramp shape trigger effect, 
based on the assumption that a gentler lee slope would allow more time for the drag to dissipate 
energy from system. Several runs, not shown here, were made with reversed topography but, at 
least for the values of Cp used herein, there was nothing to support this conjecture. 


6.4.2 Inertially unbalanced flow — (Type II) 


Further experiments were performed, with a ‘Type II’ flow balanced only against drag but sub- 
ject to topographic and inertial influence. A novel feature of these experiments was that the ini- 
tial depth was set below the height of the ranges (that is, dy) < H) and this was done to simulate 
a sea-breeze regime (or possibly a nocturnal inversion layer, or a shallow post-frontal layer). 
Layers such as these would be continuously stratified and therefore it is appropriate to deter- 
mine g’ from Eqn 14 as in the Q1'/2 model. 


In this section, we use as ‘reference’ conditions a southerly sea-breeze of 8 m/s (vp = -8 nv/s, 
Uy = 0), depth 300 m, g’ of 0.1 m/s? (roughly equivalent to a sea-breeze-layer 3 to 4°C cooler than 
the adiabatic environment) and Cp of 0.002 and E = 0. The flow is initially balanced by determin- 
ing P, and PR. from Eqn 82. The first run, R50, yields the extremely interesting result that the Cori- 
olis force gradually ‘slops’ the cooler air over the hill in the form of a supercritical density 
current. Upstream, the layer tends to be roughly horizontal at an altitude of 600 m; a value not 
greatly in excess of H itself. 
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Fig.46 Flow results for the reference conditions except that g’ = 0.1 
m/s”. Note that the supercritical flow is more extensive than 
before, although maximum speeds are decreased. 
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Fig. 47 Flow results for the reference conditions except that g’ = 0.3 m/s”. Note that the 


transitional flow is much less extensive, although maximum speeds are in- 
creased. This is equivalent to very strong gully-winds about the foothills with 
much weaker winds a few kilometres downstream. 
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Fig. 48 Flow results for the reference conditions except that Cp over the hills 
region is increased to 0.01 (though Cp over the plains is held at 
0.002). Note that the leeward jump is held back. 
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Fig. 49 Flow results for the reference conditions except that Cp over the hills re- 
gion is increased to 0.005 (while Cp over the plains is held at 0.001). 
Note that the leeward jump is held back, but less so than in Fig. 48. 
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Fig.50 Type IJ flow with initial conditions of a southerly of 8 m/s, do 
= 400 m, g’ = 0.1 m/s? and Cp = 0.002. Note that this type of 
flow is associated with a sea-breeze from the Coorong area. 


With vo = -4 m/s, there is a weaker and shallower gully-wind which occurs several hours 
later (Fig. 51). Thus, a weak sea-breeze leads to a weak or non-existent gully-wind. Obviously 
lower and smoother ranges will induce a gully-wind more easily than higher, rougher ranges. 
Runs performed for southerly sea-breezes simultaneously on both sides of the ranges gave sim- 
ilar results, but with a jump-like feature on the lee plain where the supercritical flow ran into 
the other sea-breeze. However, these are not shown here. The software as provided does not 
model this dual sea-breeze regime without a minor change to the source code. 


As a point of conjecture, it was surmised that an interfacial drag increases the possibility 
of transitional flow. If the homogeneous upper layer were to flow from left to right (easterly) 
for synoptic reasons then interfacial drag between this layer and the lower cooler layer would 
produce upslope drag. In the Southern Hemisphere, a high-pressure ridge building to the south 
at night would be consistent with these preconditions. 


It is difficult to estimate E, the interfacial drag coefficient, though on scales different to 
those of the gully-winds and for markedly more stable air, katabatic wind studies by Manins 
and Sawford (1979) show E to be of order 107. It is stressed that the following run with a 
nonzero value of E is intended only as a basis of conjecture. For vp = -4 m/s and Cp = 0.002, it 
was shown that weak downslope winds occurred (Fig. 51). However, with interfacial drag op- 
erating in such a way that w4 = 15 m/s and E = 0.001, the spillover effect is greatly enhanced 


(see Fig. 52). This suggests that the interfacial drag may indeed be a significant factor in gully- 
wind formation. 
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Fig. 51 Type II flow under the same conditions as for Fig. 50, except 
that vo = -4 m/s. 
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that E = 0.001. 











Wind speed (m/s) 


20 
10 


20 


wind speeds (m/s) 


Component 


Chapter 6 67 








6.5 Summary 


Geostrophic (and frictionless) flow perturbed by topography behaves similarly to non-rotation- 
al flow except that the upstream jump is retarded while the downstream jump is accelerated lee- 
ward. Some veering and backing occurs due to the near-conservation of the vertical component 
of potential vorticity. With increasing friction, the leeward jump is reduced in strength and re- 
tarded, sometimes to a stationary position, and the maximum wind speed is significantly re- 
duced. Presumably significant friction could inhibit the occurrence of downslope winds. 


For inertially unbalanced flow, a range-parallel sea-breeze may turn sufficiently to spill 
over the crest in a shallow supercritical-layer. This effect is reduced and delayed by increasing 
surface drag, although the details of this have not been given herein. Not surprisingly, it is also 
enhanced by interfacial drag from a favourable upper wind, although this last finding is qualita- 
tive at best. In the Southern Hemisphere (where most of this study has focussed), the presence 
of topographic barrier with a quasi-meridional orientation tends to create a southerly jet on the 
eastern flank of the ranges. The Murray Mallee area to the east of the Mount Lofty Ranges is 
known to be prone to strong sea-breezes, so this is certainly a valid result in the South Austral- 
ian case. Maybe, it could be regarded as a small scale version of the southeast Australian south- 
erly buster. 


Simulations, using realistic preconditions pertinent to the postulated archetypes, reveal 
that the model either reproduces or partially explains the main features of the Adelaide gully- 
wind. In particular, for light geostrophic wind situations, it appears that a sea-breeze spilling 
over in a shallow layer (~100 m) fits the acoustic sounder observations and balloon pilots’ expe- 
rience outlined in section 3.4. 


The model has been incorporated as a menu-style software package for practical forecast- 
ing. Operational forecasters? have experimented with the model in a ‘what if’ mode and indicat- 
ed to the author that its results are reasonably consistent with their experience. The reader is 


invited to carry out his/her own experiments using the copy of the software package provided 
with this thesis. 





3. The author, too, has many years operational weather forecasting experience. 


CHAPTER 7 


DISCUSSION AND CONCLUSIONS 


7.1 Hydraulics versus waves 


Until recently, the hydraulic theory was confined to one or two homogeneous layers. The work 
of Smith (S85 and SS87), Baines and Guest (1988) and Baines and Granek (1990) in both tank 
and numerical experiments has extended the applicability of the theory from homogeneous flu- 
ids to continuously stratified fluids. The steady-state model developed in Chapter 4 is in line 
with this approach. Recent literature (see sections 2.6 and 2.7) has expressed a preference for 
the hydraulic theory to explain the essential physics and nonlinear behaviour of downslope 
wind events. Notable exceptions are the work of Peltier and Clark (1977 and several later pa- 
pers) (PC) and of Laprise and Peltier (1989a,b,c) (LP). LP present much numerical evidence in 
support of the nonlinear wave-amplification mechanism of PC. 


For the Adelaide gully-wind and the bora, at least, the observed ratio of the height of de- 
coupling to the mountain height is about two — much less than the value of about five predicted 
by PC. It seems that such low-level decoupling is therefore due to other reasons, e.g. subsid- 
ence inversion or a shear layer. This does not rule out the PC mechanism occurring for certain 
types of flow, e.g. possibly the chinook. For the Adelaide gully-wind there is evidence (al- 
though it is not documented) that much of the flow is on occasion channelled through the gul- 
lies and valleys. One of the author’s professional colleagues has noted that gully-winds can 
often be heard flowing through a valley below his residence. Arakawa (1969) and Pettre (1982) 
have successfully used the concept of hydraulic flow through a valley to model downslope 
winds. It is difficult to see how PC’s wave-amplification theory could easily accommodate the 
channelling effect. 


The hydraulic model accounts for the postulated Type II gully-wind, where the primary 
mechanism is the Coriolis force ‘slopping’ a shallow dense layer up and over the ranges, possi- 
bly aided by drag from aloft. 


7.2 Utility of the steady-state (Q1!/,) model 


From Eqn 27 it is seen that maximum wind speeds of two to three times the upstream value are 
anticipated over the lee plains — a realistic result. Qualitatively, the model compares well with 
each of Smith (1985), Durran and Klemp (1987), Bacmeister and Pierrehumbert (1988) and PC 
and LPa,b,c. Detailed comparisons with DK in particular were very good. For instance, in the 
presence of an environmental critical-layer, both models were quite similar in predicting transi- 
tional flow. Drag predictions from Q1'/. when compared to DK and BP were extremely good. 
In contrast, drag as calculated by S85 was about double the expected value. The splitting- 
streamline height of the steady-state model tended to be 10 to 15% higher than that of S85. 
Both QI'/, and S85 predict multiple higher altitude values of dg; a result validated by DK and 
BP. The disadvantages of both these hydraulic models are that No and ug must be idealised as 
constant with height and that the existence of a critical mountain for which the transitional flow 
is self-initiated cannot be directly inferred. Baines and Granek’s quantitative approach obvious- 
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ly has many of the advantages of both models and can be used with realistic upstream atmos- 
pheres and topography. Data from tank experiments by Rottman and Smith (1989) and 
numerical experiments by DK, BP and PC show that Q1'/> is equally good at predicting do 
(Fig. 33). 


An assumption in the model is that the hydraulically-active layer may be approximated 
by an equivalent homogeneous model and that the value of g’ used for the equivalent layer is 
4| Nod f- The overall utility of the model indicates the validity of this assumption. 


The inclusion of some topographic effects, such as lower leeside plains and channelling 
between the major hills of a range, showed that valley mouths, especially those at relatively 
low elevation, are prone to chronic and severe downslope winds. This result is similar to Ar- 
akawa’s (1969) observations of the hiroto-kaze in Japan and with many observations concern- 
ing the so-called extreme katabatics off the Antarctic and Greenland coasts. 


With the remarkable exception of these extreme katabatics, a feature of all downslope 
winds is their gustiness in speed yet constancy of direction. Surge periods of 3 to 4 minutes 
have been reported for the bora, and lidar measurements of a chinook windstorm by Nieman et 
al. (1988) showed a tendency for periodicity of 4 minutes and 14 minutes. Resonant tuning 
was hypothesised between vertically and horizontally propagating waves (both of which have 
their oscillations in the vertical). Under this hypothesis directional fluctuations would not be 
anticipated — nor do they occur. Using the Q1'/, steady-state model the consequent speed fluc- 
tuations were found to be of the order of 3 minutes and 15 minutes for large-scale downslope 
winds — both values being well supported by the observations. The internal waves could be ini- 
tiated by topographic irregularities. Over ice, the lack of irregularities to initiate such waves 
could well account for the well-known constancy in the speed of the Antarctic katabatics. 


7.3 Utility of the time-dependent (QT) model 


Although this model is based on the shallow-water equations where the layers are usually as- 
sumed to be homogeneous, the results of the Q'/, steady-state model and the good agreement 
with DK’s results suggest an extension to the case of a continuously-stratified shallow layer. 
Thus, it may be used in model atmospheres where No is constant with height. 


The shallow-water equations incorporating f, CD and E were solved by a two-step Lax- 
Wendroff numerical scheme. The results are consistent with and anticipate many of the fea- 
tures of downslope winds in general and the Adelaide gully-wind in particular. 


The primary output values of the model such as maximum wind speed, zone of strong 
winds and layer depth (or inversion height) as presented are realistic. 


In simulating typical synoptic conditions expected for Type I Adelaide gully-winds (Figs 
41 to 49), the resulting leeside inversion turns out to be at about 600 m with maximum winds 
in excess of 20 m/s and these results are consistent with the data presented in section 3.2. The 
appropriate comparison for these modelled maximum winds is the observed gust speed. This is 
because, although friction has been accounted for in a parameterised manner, the actual sur- 
face wind will be significantly less than the layer-mean. Air transferred from the mid-layer re- 
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gion to produce a gust at the surface will have a speed representative of the layer-mean. In 
section 3.3, a shallow post-frontal layer crossing the ranges (in synoptic balance) evolved into a 
gully-wind with downstream depth of ~550 m, leeside depth of ~300 m and upstream depth of 
~600 m. A simulation (not shown) of this event also gave comparable values. In section 3.4 
data from an acoustic sounder and private communication from an experienced balloon pilot 
both showed that the inversion could be as low as 60 to 200 m. Such small depths and the gen- 
erally lower wind speeds are probably not explicable in terms of Type I behaviour; nor does it 
seem likely that the mechanism of PC is relevant here. They are, however, explicable in terms 
of Type II behaviour, which involves the inertial turning of a very shallow upstream layer 
(dy <H). The Type II preconditions would apply to streams such as a sea-breeze-layer or a 
shallow post-frontal layer. 


Figures 50 and 51 for Type II flow depict the evolution of this spillover due to Coriolis 
turning. With a weaker sea-breeze or greater surface drag, the spillover effect is delayed or pre- 
vented. In more speculative vein, however, it does seem that the effect of interfacial drag from 
a fast-moving upper layer would be adequate to help pull a shallow stable layer up hill even in 
the presence of surface friction. In all of these cases the subsequent leeside layer is of the order 
of 100 m deep, which is significantly shallower than those of Type I and yet quite consistent 
with the available observations. 


The sensitivity of the occurrence of transitional flow and of the jump position to upstream 
velocity, inversion height or environmental critical-layer height, stability of the layer or 
strength of the inversion, mountain height, external pressure-gradient forces and surface drag 
all account for the marked spatial and temporal variability of the Adelaide gully-wind. The in- 
clusion of non-zero Cp had the effect of bringing the jump back towards the crest. With Cp 
maximised during the day and minimised during the night, it is be expected that the Adelaide 
gully-wind, and downslope winds generally, would have a pronounced preference for occurring 
nocturnally. With Cp relatively high over the ranges, the gully-wind was inhibited. This finding 
is also consistent with the occurrence of stronger gully-winds in areas where the ranges are 
smoother and sparsely timbered. 


For Type II flow, an intriguing result is the effect of variations in Cp and H, the mountain 
height. It is apparent that increasing Cp and/or H lessens the likelihood of the spillover, so that 
rougher and/or higher parts of the ranges could be untouched by the shallow layer whilst lower 
and/or smoother parts of the ranges are subject to a shallow layer moving cross-range. In such a 
case, it follows that there would be no downslope winds over the leeslopes of Mount Lofty it- 
self, but downslope winds could be occurring on leeslopes further to the north and south. This 
may have been the situation for the nights of 2/3 and 3/4 February 1985 and 25/26 February 
1991. 


A pragmatic feature of the model is its use as a semi-interactive forecasting tool. The 
main input parameters of do, ug and g’ may be estimated from synoptic data by experienced 
forecasters. Informal feedback from the author’s colleagues indicates that the model outputs 
tend to correspond to their subjective expectations. A major problem is the correct specification 
of the initial or undisturbed conditions, since there are no upper-air data acquired routinely 
within 200-300 km upwind of the ranges. Therefore the model is difficult to either test or use in 
a thorough and objective fashion. 


CHAPTER 8 


FINAL REMARKS AND FUTURE DIRECTIONS 


A simple steady-state hydraulic model of downslope winds was developed, presented and 
found to account for many of the observed features of downslope winds. As an extension of 
this steady-state model, a time-dependent quasi layer-and-a-half model was constructed, pre- 
sented and also shown to account for several intriguing aspects of downslope winds and of the 
Adelaide gully-wind in particular. 


Several unique postulates (concerning the steady-state model as presented, mountain drag, 
turbulence intensity and jump strength as related functions of the non-dimensional mountain 
height, a resonant tuning of gravity waves to account for the surging characteristic of downs- 
lope winds, a modelled valley effect, the evolution of a downslope wind due to inertial forcing, 
the way in which drag inhibits downslope winds during daylight, and the signature of a hydrau- 
lic jump ‘written in the wind’) were proposed and framed within current knowledge and theory 
in the literature. Also robust, interactive and operationally useful software was devised. 


Not surprisingly there are still many fascinating questions to ask... and to answer. From an 
observational perspective, just how extensive are, for example, the Adelaide gully-winds, how 
far to the lee do they blow, what is their spatial variation and their association with surface fea- 
tures and conditions? The advent of automatic weather stations, loggers, computer power and 
techniques together with sophisticated and robust flying instrument packages (the aircraft oper- 
ated by the Flinders Institute of Atmospheric and Marine Science) promises much in this area. 


Noting the success of a tank model for the famous Spillane eddy in Melbourne, it is an ob- 
vious step to construct one for the Mount Lofty Ranges, especially if the results were to be com- 
pared with the data of the observational program envisaged above. 


In hindsight, it is apparent that the direction shown by Baines and Granek (1989) has the 
potential to be the simplest possible, yet realistic, model of downslope winds. Its extension to 
time-dependency will almost certainly repay the effort. Meanwhile, there are many sophisticat- 
ed mesomodel packages that will soon be easy to use on typical PCs and work stations (over 
the next 3 to 5 years). Prior to such numerical models being available to operational forecast- 
ers, it should be possible to extend the software provided with this thesis in two useful ways 
(but probably not simultaneously). Firstly, a time-dependent multi-level two-dimensional model 
would be relatively simple to construct. It would take longer to run of course (of the order of an 
hour rather than a few minutes), but would probably still be useful operationally. A second ex- 
tension would be to an x,y plane, whilst retaining the 1'/-layer vertical structure (or further ex- 
tended perhaps with 2!/2, or possibly 3'/2 layers). 


All three programs are more general than just downslope wind studies. A host of ques- 


tions relating to sea-breezes, early morning wind on the Adelaide plains, channelling and fun- 
nelling, fog advection and local eddies could be pursued. 
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APPENDIX A 


Derivation of transitional flow conditions 


As in Fig. Al, consider two-dimensional flow of water of depth d(x) over an obstacle of height 
h(x) and maximum height of H at x = 0. The flow is from left to right in the positive x direc- 
tion. The layer velocity is « = u(x) and is constant with height. Differentiation with respect to x 
is denoted by subscript x. Subscripts 0 and c denote upstream and crest values, respectively. 








Fig. Al | Schematic diagram of transitional flow over an obstacle. 


The continuity and momentum equations may be expressed as 
(ud), = 0 (Al) 


0 


uu,+g(d+h), 
eliminating u, between the equations of Al and using 


a) u 


a = — A2 
Fe (A2) 
then yields 
h, 
d= (A3) 
: Fr-] 
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Equation A3 shows that for h, > 0 and F < 1 then d, < 0 and hence u, > 0. That is, subcritical 
‘windward’ flow is accelerated. Similarly, supercritical leeward flow is accelerated. To have ac- 
celeration on both sides of the obstacle it is apparent that there must be a transition from F < 1 
to F > 1. Assuming that d(x) is a smooth function (in the sense that its first derivative exists and 
is continuous everywhere), we can show from the first equation of AJ that u,/u must be continu- 
ous and hence that u(x) is likely to be smooth (unless any discontinuities in u, and u which may 
occur happen to coincide and to cancel each other out, which is fairly unlikely from a physical 
point of view). If the topographic function A(x) is also smooth, it can be shown from A3 that 
F(x) must be everywhere continuous. Since the latter must take all values during its transition 
from F < 1 to F > 1, it follows then that there must exist some point at which F = 1. 


Since h, = 0 at the crest of the mountain, the condition for transitional flow is therefore 
that F.= 1. As F=1 atx =x,, application of L’H6pital’s rule to A3 will show that 


(A4) 


which is finite for the kind of topography considered. Integrating Eqns A1, we find that 


ugdy = u,d, (AS) 
Mo d d.+H 
ag" b= Ta ot 


and substituting A2 into AS gives 


2 H d, 12 
For 2(l=—) = 42) (A6) 
dy dy 
3% 
dy 
But from A2, 
us 
d.=— (A7) 
&§ 
and thus, 
d. ue xe 
dy ad’ a 


From AS then, 


u. = Ug. (AY) 
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Together, A8 and A9 yield 


d 3 
Fi = (=) . (A10) 


Substituting A10 into A6 gives 


F3-3F 421-2) = 0. (All) 
0 


For a homogeneous fluid overlaid by a much deeper one it may be shown that the above deriva- 
tion is valid except that g is to be replaced by the modified gravity under the Boussinesq ap- 
proximation (Baines, 1987). 


APPENDIX B 


Derivation of steady-state jump conditions 








Fig.B1 Schematic diagram of an idealised hydraulic jump. 


Consider steady-state two-dimensional flow of water with a stationary hydraulic jump as in 
Fig. B1. With subscripts 1 and 2 denoting pre- and post-jump conditions, respectively, then the 
continuity equation may be expressed as 


u,d, = Und, =m. (B1) 


By equating the momentum flux change across the jump to the pressure force across it, then it 
is evident that 


uid, - usd, = ; g(d; -d’) 3 (B2) 


Eliminating u, and uw between the three equations in B1 and B2 easily yields: 


d,— d, 


7 d,d, 


= 5 8(4— di). (B3) 


Discarding the trivial solution of d, = d), we find that the remaining root of B3 satifies 
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80 
d,+d 
m 2 1 
a ee BH) 
gd, 2d; 
from which it follows that 
(B5) 


g 
Again, the substitution of modified gravity, g’ for g, in the above equations may be made in the 
case of a shallow layer of homogeneous fluid overlaid by a deeper homogeneous layer. 


d 
a = lt Fes. 


APPENDIX C 


Derivation of steady-state mountain drag 


The drag force, r, exerted by the mountain to change airflow from what it would have been oth- 
erwise, is the difference in the combined momentum flux and pressure force acting on the 
changed layer (that is, below altitude, dj, along which p = po*) between upstream and down- 
stream sites (see e.g. S85; and indeed, the derivation given here closely follows S85, including 
nomenclature). 





| 
UNDISTURBED FLOW 


i neal" aca aaa Po 
) 


0 WELL MIXED 


REGION 











Fig.C1 Schematic diagram of homogeneous fluid crossing a two-dimensional obstacle. 


As in Fig. Cl, consider a deep homogenous layer of density p, overlaying a relatively shallow 
homogenous layer of density, p. The term PF (horizontal pressure force per unit width normal 
to the flow) at upstream site 0 in the active layer is 


PF) = pitdy= pads. (Cl) 
Similarly, at a lee side site 

PF, =pi*d, + sped (C2) 
where (from the hydrostatic equation) we have 


P; =pe +p.g (dy —d,) 4 (C3) 
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The momentum flux, MF, at each site is 


MF, = puod, (C4) 


MF, 


2 
pujd, 
Finally, the drag, r, is calculated from 
r= (PF,-PF,-T) + (MF,-MF,) (C5) 
where T denotes the pressure force on the layer from the mixed region, and is given by 
1 
T =p (dy- 41) +508 dy- dy)? : (C6) 


Now remembering that 








~. 4 
s = d, = Uo (C7) 
and that 
P.—P 
: (C8) 
p 
_ o=@. 
8 0, 
then, after some manipulation, C5 reduces to 
1 e 
r= 58(P-P,) (do—dj) — pugdy (3-1) (C9) 
1 , i 
r=38 p (do — dt) — pugdy (3-1). 
Now if we define non-dimensionalised drag as 
, _ No 
ior (C10) 
Puy 
it follows that 
,_ 243 1 a 
TF = —do( ls) +ay(S—T) (C11) 
wT Ky 


represents the non-dimensional drag for continuously stratified flow over an obstacle. Since 3 
and do are functions of A, then so too is *. 


PUBLICATIONS 


GRACE, W.J. 1991. Hydraulic jump in a fog bank. Aust. Met. Mag. 44, 205-9. 


GRACE, W.J. 1992. Modelling the downslope winds of Adelaide, South Australia. CTAC 91 
(Conf. proc. Computational Techniques and Applications 1991) (in press). 


GRACE, W.J. and HOLTON, I.C. 1990. Hydraulic jump signatures associated with Adelaide 
downslope winds. Aust. Met. Mag., 38, 43-52. 


SOFTWARE 


Files 
1 GULLY.EXE 
2 EGAVGA.BGI 


Place both files into a common directory, A:, say. Then type ‘GULLY’ and enter. The program 
is then reasonably self-evident. Use the ESC key to obtain a menu. 


All three publications and a disc containing the software file are in the back-cover jacket.* 


* Applies to original set of theses copies only. Software and publications are available from the author 
c/- Bureau of Meteorology, PO Box 421, SA, 5071 or e-mail w.grace @bom.gov.au 
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